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EXECUTIVE  SUMMARY 


This  report  is  addressed  to  analysts  with  some  background  in  probability, 
statistics,  and  stochastic  networks.  The  literature  on  activity  networks  has 
identified  shortcomings  of  deterministic  methods  such  as  the  Program  Evaluation 
and  Review  Technique  (PERT).  Much  of  the  estimation  error  of  PERT  is  due  to 
the  failure  to  adequately  treat  subnetworks  of  parallel  activities.  Statistics 
such  as  mean,  standard  deviation,  and  quantiles  of  completion  time  of  these 
subnetworks  can  be  accurately  calculated  by  the  methods  of  this  report.  There¬ 
fore,  the  errors  in  these  statistics  produced  by  deterministic  methods  can  be 
evaluated  for  specific  examples.  However,  in  the  author's  opinion,  in  most 
instances  deterministic  networking  methods  should  be  abandoned  in  favor  of 
accurate  and  comprehensive  stochastic  techniques  such  as  the  Venture  Evaluation 
and  Review  Technique  (VERT). 

The  numerical  procedures  developed  here  have  considerable  use  apart  from 
application  to  activity  networks.  The  problem  of  finding  descriptive  statis¬ 
tics  for  the  maximum  of  a  set  of  positive  continuous  random  variables  is  found 
in  the  areas  of  analysis  of  engineering  tolerances  and  of  reliability  and  main¬ 
tainability. 

Motivated  by  networking  problems,  various  parametric  analyses  are  performed. 
Goals  are  to  determine  the  sensitivity  of  the  subject  statistics  to  network 
parameters  and  to  draw  pertinent  inferences  for  activity  networks.  Some  para¬ 
meters  of  interest  are  the  number  of  random  variables  (RV's)  in  the  set,  the 
relative  size  of  the  mean  value  of  each  of  the  RV's  in  the  set,  and  the  func¬ 
tional  form  of  and  shape  parameters  for  the  probability  distributions  of  these 
RV's.  Most  of  the  examples  use  probability  distributions  having  the  domain 
zero  to  infinity.  To  meet  the  objection  that  the  domain  is  always  finite,  an 
analysis  is  made  of  the  effect  of  truncating  the  above  distributions.  A  brief 
analysis  is  also  done  on  the  effect  of  network  logic. 
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MEMORANDUM  REPORT 

SUBJECT:  Statistics  for  the  Maximum  of  Several  Positive  Random  Variables  with 
Application  to  Networking 


1.  Reference: 

References  are  designated  by  bracketed  numbers  and  are  included  in  the  foot¬ 
notes.  The  references  are  also  listed. 

2.  Background 

Deterministic*  networking  methods  for  estimating  the  time  to  complete  a 
multi-activity  project  have  existed  for  more  than  two  decades.  Methods  of  this 
sort  such  as  PERT  and  CPM,  while  useful  in  some  respects,  can  produce  large 
errors  in  estimates  of  the  mean  value  and  upper  quantiles  of  the  probability 
distribution  of  project  completion  time.  The  error  is  due  principally  to  the 
insistence  that  a  deterministic  critical  path  pass  thru  that  activity,  in  a  set 
of  parallel  activities,  which  has  the  largest  mean  completion  time**.  The 
errors  of  deterministic  networking  methods  were  identified  very  early.  For 
example,  see  Grubbs  (1962)[1]. 


*  As  opposed  to  Monte-Carlo  methods. 

**  The  deterministic  approach  replaces  the  random  activity  times  by  their 
mean  values  to  find  a  critical  path  thru  the  network.  Using  the  PERT 
assumption,  the  variance  of  the  project  completion  time  is  just  the  sum  of 
the  activity  variances  along  the  critical  path.  To  estimate  quantiles  of 
the  distribution  of  project  completion  time,  PERT  makes  the  additional 
assumption  that  the  distribution  of  this  random  variable  is  Normal. 

[1]  Grubbs,  F.;.  "Attempts  to  Validate  Certain  PERT  Statistics  or  'Picking  on 
PERT1,"  Opus.  Res.,  Vol .  10,  pp.  912  -  915,  1962. 


An  analytic  approximation  was  devised  by  Clark  (1961)[2]  to  calculate  the 
mean  and  standard  deviation  of  project  completion  time.  The  method  is  based 
on  Normal ly-distributed  parallel  activity  times.  These  times  are  permitted  to 
be  correlated,  if  activities  share  a  common  node.  Subsequent  critiques  of  PERT 
such  as  MacCrimmon  and  Ryavec  ( 1964) [ 33  and  Greer  ( 1983 ) [4] ,  attempted  to  quan¬ 
tify  the  approximate  magnitude  and  sign  of  the  estimation  errors  of  completion 
time  statistics.  By  necessity,  these  papers  use  quite  simple  examples,  from 
which  their  inferences  are  drawn.  Unfortunately,  there  is  a  great  diversity 
in  network  structure,  and  no  typical*  network  can  be  identified. 

3.  In  at  least  one  class  of  network  problems,  there  are  multiple  --  typically, 
2  to  20  --  independent,  parallel  activities  throughout  very  large  networks. 

This  type  of  network  [5]  has  been  used  by  Department  of  Army  organizations  to 
plan  for  the  reactivation  of  inactive  ammunition  production  facilities.  Moti¬ 
vation  for  the  work  reported  here  followed  discussion  with  one  of  the  authors 
of  [5]  (Moeller)  about  the  estimation  errors  associated  with  deterministic 


[2]  Clark,  C.E  "The  Greatest  of  a  Finite  Set  of  Random  Variables,"  Opns.  Res. 
Vol.  9,  No.  2,  pp.  145  -  162,  March  -  April  1961. 

[3]  MacCrimmon,  K.R.  and  Ryavec,  C.A.  "An  Analytical  Study  of  the  PERT 
Assumptions,"  Opns.  Res.,  Vol.  12,  No.  1,  pp.  16  -  37,  January  -  February 
1964. 

[4]  Greer,  W.R.  Jr.  "Why  Doesn't  PERT  Work?,"  Resource  Mgmt.  Journal,  pp.  27  - 
31,  Summer  1983. 

*  There  are,  of  course,  classes  of  problems  which  yield  similar  networks. 
However  exceptions  are  manifest  to  generalizations  such  as  "activity  net¬ 
works  have  few  independent  parallel  activities." 

[5]  Matheiss,  T.H.,  Moeller,  G.  and  Kilar,  J.  "Improving  Industrial  Readiness 
with  Venture  Evaluation  and  Review  Technique  (VERT),"  Interfaces,  Vol.  12, 
No.  1,  pp.  21  -  26,  February  1982. 


networking  methods  applied  to  this  type  of  problem.  It  appears  that  the  defects 
of  deterministic  networking  methods,  tho  adequately  reported*,  have  not  been 
appreciated  by  all.  The  beginning  of  the  present  effort  was  an  attempt  to 
quantify  the  error  of  PERT  completion  time  for  a  subnetwork  of  multiple  paral¬ 
lel  activities. 

4.  An  approach  to  all  stochastic  networks  which  avoids  the  restrictive  assump¬ 
tions  of  PERT  is  called  VERT,  for  Venture  Evaluation  and  Review  Technique.  A 
recent  expository  paper  by  Moeller  and  Digman  ( 1981 ) [6]  describes  the  modeling 
technique  and  illustrates  this  with  an  example  of  planning  in  the  electric 
power  industry.  The  textbook  on  VERT  [7]  provides  greater  detail. 

5.  Objectives  and  Scope 

The  class  of  problem  posed  by  parallel  subnetworks  of  activities  of  random 
duration  is  the  following.  Given  a  set  of  k  positive,  continuous  random  vari¬ 
ables,  each  having  a  possibly  unique  distribution,  what  are  the  values  of 
statistics  associated  with  the  largest  member  of  this  set?  The  statistics  of 
interest  are  the  mean,  standard  deviation,  and  quantiles  of  the  distribution 
of  the  largest  value.**  Objectives  of  this  report  are:  (a)  to  develop  a 
general  procedure  for  obtaining  this  mean  and  standard  deviation  and  (b)  to 
derive  some  analytic  results  suitable  for  assessing  the  computational  error  of 


*  In  fact  some  authors  ([3]  and  [4])  have  advocated  patching  up  PERT  via 
analytic  corrections. 

[6]  Moeller,  G.L.  and  Digman,  L.A.  "Operations  Planning  with  VERT,"  Opns.  Res., 
Vol  29,  No.  4,  pp.  676  -  697,  July  -  August  1981. 

[7]  Lee,  S.M.,  Digman,  L.A.,  and  Moeller,  G.L.  Network  Analysis  for  Management 
Decisions,  Kluwer-Nijhoff  Pub.,  Hingham,  MA,  c.  1981. 

**  A  method  for  obtaining  these  statistics  for  this  problem  has  greater  applic¬ 
ability  than  just  to  activity  networks.  Other  applications  lie  in  the 
areas  of  tolerance  stackup  (where  parallel  elements  exist)  and  in  reliability 
and  maintainability  analysis. 
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the  numerical  method.  Additionally,  considering  the  original  problem  context, 
certain  numerical  generalizations  are  desired  which  can  be  applied  to  activity 
networks. 


6.  Method 


Prior  to  discussing  methods,  some  problem  nomenclature  will  be  useful.  The 
probability  density  function  (p.d.f.)  and  its  integral  the  cumulative  distribu¬ 
tion  function  (c.d.f.)  for  each  of  the  random  variables  (RV's)  in  the  set 
{x.,1  <  i  _<  k }  are  denoted,  respectively,  by  f. (x)  and  F. (x).  These  are  the 
primary  inputs  or  problem  ingredients.  It  is  assumed  that  the  mean  and  variance 
of  xi  --  E[x.]  and  V[x-]  --  are  readily  available.  The  random  variable  of 
interest  is  denoted  by  z. ,  where 


zk  =  max(xp  x^,  ...,  x^,  ...,  x^),  0£x^  <  ® 

The  p.d.f  and  c.d.f.  of  z^  are  denoted  by  gk(z)  and  G^(z).  Given  that  an  expres¬ 
sion  for  g  (z)  can  be  derived,  one  can  obtain  an  analytic  expression  for  the 


j  th  origin  moment  a  -  ( k )  by 

fCO  ^  'J 


k)  =  j  r  gk(z)dz,  0  1  j 


From  this  result  the  mean  and  variance  of  zk  are,  immediately, 

E[zk]  -  a, (k ) 

n 

VL2k]  -  a2(k)  -  aj(k) 

with  standard  deviation  of  zk  equal  to  An  expression  for  g k ( z )  is 

easily  obtained  for  certain  types  of  F-(x)  distributions.  Analytic  results  for 
examples  of  this  sort  are  obtained  by  the  above  method.  These  results  are 
derived  in  Annex  A.  Also  included  in  this  annex  are  probability  arguments 
leading  to  general  procedures  for  calculating  E[zk]  and  V[zk],  which  can  be 
easily  implemented  in  a  computer  program.  From  a  pragmatic  point  of  view,  it 
is  immaterial  whether  numerical  results  are  obtained  by  evaluating  a  closed- 
form  expression  or  by  following  another  numerical  procedure,  providing  the 
latter  is  not  computationally  too  expensive  and  yields  a  sufficiently  small 


numerical  error.  A  computer  program  was  developed  to  implement  the  procedure 
derived  in  Annex  A.  The  program  listing  is  displayed  in  Annex  B.  The  method 
presented  here  has  the  important  practical  advantage  of  being  free  of  restric¬ 
tions  on  the  form  of  F^(x).  For  simplicity,  numerical  examples  of  this  method 
were  calculated  for  some  familiar  two-parameter  distributions  --  gamma  and 
Weibull  distributions.  A  parametric  analysis  is  conducted  which  systematically 
examines  the  effects  of  the  number  (k)  of  RV's  in  the  set  and  of  the  parameters 
of  these  probability  distributions. 

7.  As  indicated,  the  primary  emphasis  in  this  report  is  on  the  maximum  of  a 
set  of  k  random  variables.  Ordinarily,  the  logic  of  an  activity  network  cor¬ 
responds  to  this  problem.  Passage  to  other  activities  downstream  of  a  set  of 
parallel  activities  is  conditional  upon  completing  all  of  the  parallel  activi¬ 
ties.  Occasionally,  network  logic  permits  passage  when  only  k  of  n  activities 
are  complete.  _If  the  probability  distribution  of  all  parallel  activities  is 
the  same,  the  latter  problem  reverts  to  a  problem  in  order  statistics.  Guenther 
( 1977 ) [9]  presents  an  easily  applied  numerical  technique  for  evaluating  the 
c.d.f.  of  the  kth  ordered  (in  algebraic  magnitude)  statistic  in  a  set  of  n. 

The  method  of  [9]  is  exploited  here  for  this  special  case.  Details  are  pre¬ 
sented  at  the  end  of  Annex  A. 

8.  Numerical  Results 

Following  the  derivation  of  equations  in  Annex  A  some  numerical  examples 
are  considered.  The  random  variables  (x^)  for  all  examples  are  scaled  so  as 
to  facilitate  comparison  between  examples.  The  parameters  in  F^ (x)  are 
selected  so  as  to  make  the  largest  over  i  of  E[x.],  1  <  i  <  k,  equal  to  unity. 

In  fact,  one  may  as  well  order  the  x^  in  order  of  decreasing  mean  value.  This 
scaling  does  not  reduce  the  generality  of  our  approach.  (One  can  always  con¬ 
vert  between  units.)  It  does,  however,  permit  one  to  compare  the  value  of 


[9]  Guenther,  W.C.  "An  Easy  Method  for  Obtaining  Percentage  Points  of  Order 
Statistics,"  Technometrics,  Vol .  19,  No.  3,  pp.  319  -  321,  August  1977. 
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E[zk]  with  the  PERT-est imated  value,  which  is  always  unity.  Altho  the  method 
for  calculating  the  statistics  of  z^  does  not  require  the  location  parameter 
of  F.(x)  to  be  zero,  all  numerical  results  presented  here  assume  that  x^  is 
bounded  from  below  by  zero.  It  is  recognized  that  this  assumption  may  be 
unrealistic  for  certain  applications. 

9.  The  first  examples  treat  the  case  in  which  F^(x)  is  exponential  with  rate 
parameter  x^.  An  analytic  solution  exists  for  this  case.  This  permits  cal¬ 
culation  of  the  computational  errors  in  E[z ^3  and  V[z^]  associated  with  the 
numerical  method.  Results  for  the  special  case  in  which  X-  is  unity  for  all 

i  are  shown  in  Table  E.l.  The  c.d.f.'s  of  z^,  with  k  as  a  parameter,  are 
shown  in  Figure  E.l  for  this  case.  Plots  of  Gk(z)  for  all  the  examples  have 
been  made  on  Normal  probability  paper.  Advantages  to  plotting  in  this  manner 
are:  (a)  departures  from  a  straight  line  (Normality)  are  evident  and  (b)  the 
values  of  the  c.d.f.  in  each  tail  can  be  accurately  plotted  and  read.  Results 
for  other  examples  involving  exponential  F-  (x)  are  shown  in  Tables  E.2  and 
E.3.  Figure  E.2  illustrates  the  behavior  of  Gk(z)  as  k  increases,  for  the  case 
in  which  X1  =  1  and  Xi  =  1.2X._lf  i  >  1.  In  this  case  each  RV  added  to  the 
set  has  a  mean  value  that  is  1/1.2  of  its  predecessor.  Convergence  of  quantiles 
in  the  upper  tail  of  Gk(z)  is  evident.  Additional  observations  concerning 
this  and  other  examples  are  found  in  Annex  A. 

10.  Accuracy  of  the  numerical  method  is  displayed  for  two  examples  in  Tables 
E.4,  E.5,  and  E.6.  Accuracy  is  a  function  of  the  step  size  used  in  numerical 
integration  and  of  the  method  of  quadrature  used.  Two  quadrature  schemes  are 
employed  for  comparison:  rectangular  and  Simpson's  rule  [11].  The  latter  is 
preferred  on  the  basis  of  computational  efficiency,  altho  both  schemes  are 
easily  implemented  and  yield  satisfactory  accuracy.  Methods  for  obtaining 

E[ z kJ  and  V[ z^]  based  on  numerical  evaulation  of  the  integral 


[11]  Bennett,  A. A.,  M’lne,  W.E.  and  Bateman,  H.  Numerical  Integration  of 
Differential  Equations,  Dover  Pub.,  New  York,  NY,  c.  1956. 
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j  zJgk(z)dz,  j  =  1,  2 

are  not  recommended  due  to  their  relatively  poor  accuracy  vis  a  vis  the  method 
of  Annex  A. 

11.  Some  parametric  analyses  were  performed  using  the  number  (k)  of  RV's  in 
the  set  as  a  parameter  and  with  F^(x)  having  two  functional  forms  --  gamma  and 
Weibull.  The  effect  of  k  on  ( z )  for  a  Weibull  distribution  with  shape  para¬ 
meter  2  is  shown  in  Figure  E.3.  The  shape  parameter  (6)  of  these  distributions 
influences  both  the  degree  of  dispersion  and  the  skewness.  With  the  constraint 
that  E[x.]  =  1,  1  1  i  1  k,  the  shape  parameter  of  the  F^x)  distribution  was 
changed  systematically  to  determine  the  effect  on  E[zk].  Results  for  the 
gamma  distribution  are  shown  in  Table  E.7.  Comparable  results  for  F^x) 

Weibull  are  shown  in  Table  E.8.  The  effects  upon  the  coefficients  of  varia¬ 
tion  and  of  skewness  of  x^  due  to  changes  in  shape  (8)  are  different  in  the 
gamma  and  Weibull  distributions.  These  differential  effects  are  illustrated 

in  Table  E.9.  Even  when  the  two  values  of  8  are  chosen  to  give  the  same  coef¬ 
ficient  of  variation  of  x^,  one  should  expect  different  values  of  E[zk]  in 
the  gamma  and  Weibull  cases.  In  fact,  such  differences  are  observed.  Stated 
differently,  a  difference  exists  between  Efz^]  when  F^(x)  is  gamma  versus 
E[zk]  when  F^x)  is  Weibull  with  the  same  coefficient  of  variation.  However, 
this  difference  is  quite  small*  whenever  the  skewness  of  F^(x)  >  1  and 
k  £  6.  Under  these  conditions  the  form  of  the  c.d.f.  of  x^  is  not  important 
when  calculating  E[z^] . 

12.  Another  numerical  study,  applicable  to  activity  networks,  deals  with 


*  A  difference  in  these  cases  is  about  ±0.02  or  less  for  E[z^]  ^  1.7  and 
SD  zk  ^  0.4.  To  detect  a  difference  of  this  magnitude  via  stochastic  simu¬ 
lation  would  require  a  Monte-Carlo  sample  greater  than  500  for  the  standard 
error  of  the  estimate  to  be  less  than  about  0.02. 


the  upper  limit  on  the  range  of  x^.  The  previous  types  of  distributions  con¬ 
sidered  for  F^(x)  were  defined  on  the  semi-infinite  domain  (Q,  <*>).  This  was 
done  for  analytic  convenience.  However,  in  practice  some  mechanism  will  act 
to  truncate  x^  from  above.  One  may  ask  what  the  effect  of  truncating  F..(x) 
at  some  large  quantile  is  on  the  mean  and  standard  deviation  of  z^.  It  is 
clear  that  distributions  of  x^  which  exhibit  large  positive  skewness  would  be 
most  sensitive  to  truncation  for  our  problem.  Therefore,  an  exponential 
c.d.f.  was  used.  This  distribution  was  truncated  at  the  0.99  and  at  the  0.999 
quantile.  The  value  of  A.  was  adjusted  so  that  E^]  is  the  same  (=1)  for  all 
cases.  Results  of  these  calculations  are  displayed  in  Table  E.10.  One  obser¬ 
vation  of  interest  is  that  the  mean  of  z^  is  much  less  affected  by  truncation 
than  the  standard  deviation  of  z^.  Truncation  at  even  the  0.99  quantile  does 
not  have  a  remarkable  effect  on  E[z^] .  For  example,  E[Zg]  is  2.45  without 
truncation;  is  2.43  when  truncation  is  at  the  0.999  quantile;  and  is  2.36  when 
truncation  occurs  at  the  0.99  quantile.  However,  the  corresponding  standard 
deviations  of  z^  are,  respectively,  1.22,  1.15,  0.97. 

13.  Conclusions 

The  deterministic  approach  to  activity  networks,  which  replaces  random 
activity  times  by  their  mean  values,  can  lead  to  large  underestimates  of  the 
mean  completion  time  for  subnetworks  of  parallel  activities.  For  specific 
subnetworks  of  parallel  activities,  the  methods  of  this  report  can  be  used  to 
estimate  the  PERT  error  in  the  mean,  standard  deviation,  and  quantile  of  the 
completion  time.  This  approach  is  preferable  to  numerical  general izations, 
since  activity  networks  are  quite  diverse.  The  PERT  assumption  regarding  the 
Normality  of  the  completion-time  probability  distribution  is  grossly  wrong  for 
(sub)networks  having  multiple  parallel  random  activities,  each  of  which  is 
positively  skewed.  In  those  instances  where  network  structure  and  logic  are 
quite  complex,  it  seems  practical  to  use  a  stochastic  networking  technique 
such  as  VERT  rather  than  attempt  to  "patchup"  PERT. 

14.  The  method  derived  to  treat  the  problem  of  networks  of  parallel  random 
activities  is  considerably  more  general  than  this  application  may  suggest.  The 


method  derived  and  used  in  this  report  is  quite  accurate  for  calculating  the 
means,  standard  deviations,  and  quantiles  of  the  maximum  of  a  set  of  k  posi¬ 
tive,  continuous  random  variables.  In  applications  involving  activity  networks, 
the  computational  errors  of  the  method  are  completely  negligible. 

15.  Certain  quantitative  generalities  have  been  induced  from  specific  numeri¬ 
cal  examples.  It  is  noted  that  when  the  mean  of  one  of  the  exponential  RV's 
(x.)  in  a  set  of  k  is  greater  than  the  others,  the  standard  deviation  of  zk 

(=  max(x.))  will  eventually  decrease  with  increasing  k.  An  implication  for 
i  1 

activity  networks  is  the  following:  As  the  number  of  parallel  activities 
increases,  the  standard  deviation  of  the  completion  time  will  actually  decrease 
beyond  a  certain  point.  This  contradicts  one  of  the  PERT  assumptions.  For 
positively  skewed  distributions  of  x^  —  such  as  gamma  and  Weibull  F^(x),  the 
mean  of  zk  for  k  <  6  is  not  very  sensitive  to  the  form  of  F^x)  provided  the 
distribution  parameters  are  chosen  to  yield  the  same  first  two  statistical 
moments.  This  fact  makes  consideration  of  the  precise  form  of  the  distribution 
of  xi  somewhat  academic  for  parallel  activity  networks. 

16.  Consider  progressively  adding  RV’s  (x^)  to  the  set;  i.e.,  increasing  k, 

in  such  a  manner  that  each  xi  has  a  smaller  mean  than  its  predecessor.  In  this 
case  the  upper  tail  of  the  c.d.f.  of  the  maximum,  Gk(z),  is  insensitive  to  k 
above  a  certain  point.  For  example,  for  E[x..]  =  1.2  E[x..+1],  no  appreciable 
change  occurs  in  the  0.9  quantile  for  k  >  3.  The  implication  of  this  result 
for  activity  networks  is  that  adding  more  activities  to  a  parallel  network  may 
not  change  the  low-risk  forecasted  completion  time. 

1', .  The  use  of  semi-infinite  c.d.f.'s  for  x^  instead  of  truncated  c.d.f.  's 
with  the  same  mean  is  justified  if:  (a)  the  focus  of  interest  is  on  E[zk] 
and  (b)  the  truncation  point  exceeds  the  0.99  quantile.  In  some  cases  network 
logic  may  require  that  only  k  of  n  (k  <  n)  parallel  activities  need  to  be  com¬ 
pleted  before  passing  this  point  in  the  network.  Even  tho  n  may  be  large,  viz. 

>  6,  the  difference  between  mean  completion  times  for  the  two  cases  n  of  n 
versus  k  of  n  can  be  remarkable.  This  fact  suggests  that  particular  attention 
be  paid  to  characterizing  network  logic  for  parallel  activities.  Implementation 
of  diverse  logical  forms  is  facilitated  by  using  stochastic  networking  methods 
such  as  Venture  Evaluation  and  Review  Technique  (VERT). 
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ANNEX  A 


Mean  and  Standard  Deviation  for 
the  Maximum  of  Several  Positive 
Continuous  Random  Variables 


This  annex  derives  formulas  for  calculating  the  statistical  moments  of  the 
largest  random  variable  in  a  set  of  k  positive  random  variables  (x^  1  £  i  £  k). 
Application  is  made  to  an  example  in  which  x^ ,  1  £  i  £  k,  are  exponential  random 
variables  (RV's)  from  distributions  having  different  rate  parameters.  This 
example  is  sufficiently  tractable  to  allow  a  closed-form  solution  for  the 
statistical  moments.  These  exact  values  are  used  to  evaluate  the  accuracy  of 
numerical  procedures,  which  are  suitable  for  a  more  general  case.  Specific 
examples  are  displayed,  and  some  general  inferences  are  made  from  them.  The 
examples  are  chosen  for  their  applicability  to  networks  of  parallel  activities 
which  must  all  be  completed  for  passage  thru  the  network.  It  is  noted  that 
this  type  of  problem  is  a  special  case  of  the  problem  in  which  passage  thru 
the  network  requires  the  completion  of  k  of  n,  k  £  n,  activities. 

To  start,  consider  two  positive  continuous  random  variables  Xj  and  x2 
having  probability  density  functions  (p.d.f.'s)  f^(x)  and  f2(x),  respectively. 
The  associated  cumulative  distribution  functions  (c.d.f.'s)  are  denoted  F^x) 
and  F2(x),  with  a  domain  of  x:  (0  £  x  <  °°). 

Define 

y  =  max  (x^,x2) 

with  p.d.f .  f  (y)  and  c.d.f.  Fy(y).  The  one's  complement  of  F.  is  denoted 

T.  =  1  -  F.  ,  i  *  1,2  .  (2) 

Invoking  the  definition  of  y  in  (1)  and  probability  arguments,  one  can  state 


f  (y)  »  F1(y)f2(y)  +  F2(y)fj(y) 


A-l 


Suppressing  functional  notation,  this  expression  can  be  written  as 
-  d(FjF2) 


Fy  “  .  (3c 

Using  (2)  with  (3a), 

fy  =  fl  +  f2  "  ^lf2  +  ‘  (4 

The  last  expression  can  be  used  to  find  a  simple  recursive  equation  for  the 
expectation  of  y: 

E[y]  =  f0  yFy(y)dy  .  (5 

From  (4)  and  (5), 

fOO  (OO  »QO  *00 

E[y]  =  J0  yfx(y)dy  +  J0  yf2(y)dy  -  Jo  yF1(y)f2(y)dy  -  J0  yF2(y)f1(y)dy  . 
Recalling  the  definitions  of  f1  and  f^, 

E[y]  =  Etxj]  +  E[x2]  -  \0  yF1f2dy  -  \Q  yF^dy  .  (6b 

Unng  integration  by  parts, 

*QO  -CO  .00 

-  J0  yF1f2dy  =  Jo  yF2fxdy  +  JQ  F^dy  -  E[Xj]  . 

Combining  this  result  with  (6b)  produces 

E[y]  =  E[x2]  +  J0  F1F2dy  .  (7a 

One  can  obtain  an  alternative  expression  by  using  symmetry  arguments  and  by 
exchanging  indices,  or  by  the  following  argument.  From  (7a)  using  complements 
of  Fj  and  F2, 

E[y]  =  E[x2]  +  I0  ( F ^  -  F2  +  F1F2)dy 


E[y]  =  EfXj]  +  \Q  F^dy 

The  second  moment  of  y  with  respect  to  the  origin, 


can  be  obtained  from  (4): 

E[y2]  =  E[x2]  +  E[x|]  -  {o  y2F1f2dy  -  y^T^dy  .  ( 

This  expression  is  analogous  to  (6b)  for  the  first  statistical  moment.  Inte¬ 
grating  the  first  integral  in  (8)  by  parts  and  combining  terms  gives 

f°° 

E[y2]  =  E[x2]  +  2Jo  yF^dy  (9 

or 

E[y2]  =  E[x2]  +  2.^  yF^dy  .  (9 

Exampl  e 

For  a  specific  example  of  the  above  theory,  suppose  that 
Fj(x)  =  1  -  e_XlX  (E. 

and 

F2(x)  =  1  -  .  (E. 

From  (7a),  the  expected  value  of  y  is 
E[y]  =  A21  +  |“  e"Aiy(l  -  e‘X2y)dy 
or 

E[y]  =  Ai*  +  A2^  -  (Ai  +  A2 )  ^  .  (E. 

We  will  return  to  this  example  for  extension  later. 

The  relation  between  distribution  functions  in  (3)  for  the  maximum  of  two 
RV's  can  be  generalized  to  the  max  of  k,  as  follows.  Define  a  positive  RV  xk 
with  p.d.f.  and  c.d.f.  denoted  by  f^  and  F^.  Also,  define  the  RV  z^: 

-  max(xj ,  x2,  •«.,  x,j,  x^)  *  (1 

with  p.d.f.  and  c.d.f.  denoted  by  gk  and  G^,  respectively. 

Note  that 


Then,  (3)  yields 

Wz)  ■  V2'WZ)  *  Fk+i(z)9k(z) 


(12a) 


Wz>  -  Gk(z)Fk+i(z) 


(12b) 


Gk+i‘z>  ■  nj=!FJ<z> 


(12c) 


Notice  that  both  the  p.d.f.  and  c.d.f.  of  zk  can  be  obtained  recursively. 
Further,  there  is  no  requirement  that  all  the  F.  be  identical,  as  in  the  case 

J 

with  order  statistics. 

In  (12)  plays  the  role  of  F^  and  F^+1  plays  the  role  of  F 2  in  (3a).  A 
similar  exchange  of  variables  in  (7a)  produces  the  following  relation  for  mean 
values  of  zk 

E^zk+1^  =  E^xk+1^  +  °  ^13a^ 

If  "Fj<+j  has  a  simpler  form  than  G^,  the  following  analog  of  (7b)  may  be 
preferred 

E[zw]  ■  E[zk]  ♦  f0  Gk(y)Fk+1(y)dy  .  (13b) 

A  generalization  of  (9b)  for  the  second  origin  moment  is 


E[<+1]  -  E[zJ 


h  +  2j0  yGk(y)Fk+1(y)dy,  k  >  1. 


The  variance  of  zk+^  can  be  obtained  from  the  first  and  second  origin  moments 
via: 

«zk+1]  =  E[zk+lJ  -  (ECzk+l3>2  .  05) 

Note  that  (13),  (14),  and  (15)  provide  the  basis  of  a  numerical  procedure  for 
calculating  the  mean  and  variance  of  zk  recursively.  Computational  experience 
with  this  procedure  indicates  that  somewhat  better  accuracy  is  obtained  than 


A-4 


is  possible,  at  the  same  integration  step  size*,  with  a  straightforward  inte¬ 
gration  of  G.  via 

,  K  . .. 


E[zk]  - 

J0  Gk(y)dy 

(16) 

■ 

2^  yG(y)dy  . 

(17) 

These  general  results  can  be  used  to  extend  the  example,  started  following 
(9).  Suppose  that  a  random  variable  x3  also  has  an  exponential  distribution: 

f3(x)  =  A3e-X* 
and 

F3(x)  =  e"X3X  .  (E.4) 

From  (12),  with  G^  =  F^  , 

g2(z)  -  (1  -  e‘x,z)X2e-X22  +  (1  -  e'X!2)Aje-x'z 
or 

g2(z)  =  Aie'XlZ  +  A2e"A2Z  -  (Aj  +  A2)e"(Xl+X2)z  .  (E.5) 

Whence, 

G2(z)  =  1  -  e"XlZ  -  e"XzZ  +  e"(Xl+Xz)z  .  (E.6) 

Since  E[y]  in  (E.3)  is  the  same  as  E[z2],  in  this  notation,  (13b)  gives 

E[z3]  =  AJ1  +  A'1  -  (Aj  +  A2)_1  +  fQ  e'X3y(l  -  e'Xiy  -  e"Xzy  + 
e-(A1+A2)y)(jy 


*  The  same  type  of  quadrature  formula  is  also  assumed.  For  computational 
efficiency  Simpson's  rule  [11,  p.  31]  is  recommended. 

[11]  Bennett,  A. A.,  Milne,  W.E.,  and  Bateman,  H.  Numerical  Integration  of 
Differential  Equations,  Dover  Pub.,  New  York,  NY,  c.  1956. 
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VIVV  I1  >-.■ j'  j;  J'V  -■ 


or 


E[z,]  =  AJ1  +  A’1  +  A'1  -  +  \z)’1  -  (A,  +  A,)'1  -  (A,  +  A,)'1  + 


(Ai  +  A^  +  A3) 


1 


-1 


(E.7) 


The  p.d.f.  of  z3>  93(2).  is  obtained  from  (12)  using  (E.4,  E.5,  E.6). 
g3U)  =  (1  -  e-A|2  -  e“X2Z  +  e-(A,tXi)z)X3e-Xlz  ♦  (1  -  e’*32) 

(X1e'X|Z  +  X2e'X2Z  -  (Xj  t  X2)e'(XltX!  )z)  . 

Simplifying, 

g3(z)  =  Aie"XlZ  +  A2e“X2Z  +  A3e_A3Z  -  (Aj  +  A2)e"(Al+A2)z  - 

(Aj  +  A3)e‘(Al+A3)z  -  (A2  +  A3)e"(A2+A3)  +  (Aj  +  A?  +  A3) 

-( Ai+A2+A3)z 


or 


g3(z>  >  $  Aj  exp  -  (Xj2)  -  k;j  jfj  (Xj  *  V  exp  -  (Xj  +  A^z  ♦ 

j5f  xj  exp  -  (/i  xj)z  • 


(E.8) 


The  expression  for  G3(z)  is  simplified  if  one  uses  the  following  notational 

convention.  Let  E(r,z)  denote 

i  -rz 
1  -  e 

Then,  it  is  seen  that  each  of  the  terms  in  (E.8)  is  an  exponential  density 
so  that 

G3<z>  -  jfil  «xj.  *>  -  kSj  EOj  *  V  *>  *  Efjij  Aj.  z)  .  (E.9) 

Notice  that  when  the  x^  random  variables  in  (10)  are  all  exponential,  the 
density  is  composed  of  sums  of  exponential  densities.  This  is  seen  in 


(E.5)  and  (E.8).  In  the  recursive  relation  for  g^,  equation  (12),  one  notes 


for  the  specific  example  being  considered  that 

'\+lz  "Xk+12 

9k+l^  =  9k^  +  Gk^z^k+le  "  e  9k^ 


(E ,10) 


The  last  two  terms  in  (E.10)  are  seen  to  contribute  additional  exponential 
density  terms  to  those  of  (z ) .  Specifically  for  k  =  3,  using  (E.9), 

94(Z)  -  g3u)  *  v^iyi  -  y'x"2k?]  +  + 

-  e'A'2jS?xj  “p  ■  <Xjz>  * 

*‘x'2k^  jii(xj *  xk>  «p  -  <xj  *  V*  ■  *‘x‘2ji?xj  «p  - 


After  some  manipulation  this  expression  simplifies  to 

94(z)  =  j5ixj  “P  -  <y>  -  k5j  jst<xj  *  V  exP  ■  <xj  +  Vz  + 

k*l  jSi  ilJ'x>  +  xj  *  xk>  w  -  <xi  *  xj  +  Vz  - 

«p  -  lj*ixj)2  • 

Using  the  notation  of  equation  (E.9),  the  c.d.f.  of  z4  can  be  written 


(E.ll) 


(E.12) 


g4(z)  -  -  kjJ  *  \.z)  *  kSj  ,5?  ,I?E(N  *  *  y*>  - 

E(d£ixJ*z)  •  (E,13) 

Since  G.(z)  is  a  sum  of  exponential  probabilities,  where  the  typical  term 
K  -1 

E(r,z)  has  an  expected  value  of  r  ,  one  can  immediately  write 


Efz.1  =  A:1  -  E*  .eJ(A.  +  \.)~l  +  .£*.  .11  .E?(X.  +  X.  +  A  ) 
L  4J  j=l  j  k>j  j=lv  j  k'  k*j  i=lv  i  j  k' 


4  v4  r4/ 


(.^X.f1 
VJ=1  r 


(E.14) 


The  expected  values  of  z^,  z^,  and  of  z^  --  in,  respectively,  (E.3),  (E.7). 
and  (E.14)  --  in  this  example  all  conform  to  a  common  formulation:  Combine 
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'MU  ■*! 


^  V 


£ 

i  .  • 

K- 

i  . 


the  rate  parameters  associated  with  the  k  exponential  random  variables  by 

summing  sets  of  terms,  with  sets  alternating  in  sign.  The  first  set  is  just 

the  sum  of  all  k  inverse  rates.  The  second  set  is  negative.  This  second 

set  is  the  sum  of  the  inverses  of  the  sum  of  unique  pairs  of  parameters.  The 

third  set  --  if  it  exists  —  is  the  sum  of  the  inverses  of  the  sum  of  unique 

triples;  i.e.,  combinations  of  k  values  of  A^  taken  three  at  a  time.  The 

last  set  consists  of  a  single  term  --  the  inverse  of  the  single  combination 

of  k  values  of  A.  taken  k  at  a  time.  With  this  formulation  the  expectation 

of  zc  is  immediately  written  as 
t> 

5  A  ~  1  r»5  r*  5  /  A  A  \  — 1  .  r>  5  V  5  \  »  A 


e[z5]  -  iJft1  -  fii  ♦  V'1  *  k?j  iVi  fSi<N  +  xi  *  V 
A  k?j5  A  iii'N  V'1 +  <i5iV‘ 

The  c.d.f.  of  zg  is,  by  induction. 


( E .  1 5 ) 


g5(z)  -  i£fE(Ai.z)  -  $  iifEC^i  +  >j.i)  -  k5j  j5i  +  xj  +  Vz> 


l5j  M  i=lE(Xi  +  xj  +  Xk  +  Ai»z)  +  E(ii1xi.z) 


(E. 16) 


The  second  statistical  moment  of  G k(z)  with  respect  to  the  origin,  E[z2], 
on  be  obtained  directly  from  Gfc(z)  by  using  the  fact  that  the  second  origin 
moment  of  E(r,z)  is  2r-2.  Then  from  E[z2]  and  E[zk],  the  variance  of  zk  is 
written,  from  (15),  as 

V[zk]  -  £[z2]  -  (E[zk])?  • 

Thus,  from  (E.5), 

.-2  ny~2.  o/^  .  1  \~2 


ECz|]  =  +  2A'4  -  2(Aj  +  A2) 

From  (E.  3)  and  (E.17)  and  using  (15), 

V[z2]  =  A'2  +  A"2  -  3(Ax  +  A2)-2 


From  (E.9), 
2- 


E[z3]  '  2j>lAj2  '  2k5j  +  V"‘  *  2(jIlV 


-2 


i-2 


(E.17) 


(E.18) 


(E.19) 
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One  observation  of  interest  from  Table  E.l  is  that  the  coefficient  of  variation 
of  zk  decreases  as  k  increases.  In  this  example  the  mean  increases  more  rapidly 
with  k  than  does  the  standard  deviation. 


To  demonstrate  the  effect  of  differences  in  the  rate  parameters 
(i  =  1,  . ..j,  k),  consider  the  following  numerical  example.  Let  ^  =  1,  and 
for  all  succeeding  values  of  let  .  The  values  of  ECz^]  and 

E[z^]  must  be  calculated  via  the  formulas  preceding  (E.22).  Results  are 
tabulated  in  Table  E.2. 


Coef  Var  z. 


L.0000 

1.0000 

L.3788 

1.0366 

L . 5593 

1.0182 

L.6514 

0.9948 

L.6992 

0.9763 

L .7237 

0.9637 

1.7359 

0.9560 

Note  that  E[zk]  appears  to  be  near  an  asymptote  for  k  =  7.  In  this  instance  the 
coefficient  of  variation  diminishes  even  more  rapidly  with  k  than  it  does  in  the 


case  in  which  all  ^  =  1. 


In  contrast  to  the  above  examples,  consider  a  case  in  which  all  the  random 
variables  (x^)  in  the  set,  save  one,  have  the  same  mean.  The  exceptional  RV 
has  a  greater  mean.  A  numerical  example  of  this  case  is  shown  in  Table  E.3. 


TABLE  E.3 


STATISTICAL  MOMENTS  OF  THE  MAXIMUM  OF  A  SET 
OF  EXPONENTIAL  RANDOM  VARIABLES  x.,  1  <  i  <  k, 
WHERE  E[xt]  =  1  AND  Efx^  =  0.79,  i  >  1 


k 

E[zk] 

Std  Dev  z k 

Coef  Var  z^ 

1 

1.0000 

1.0000 

1.0000 

2 

1.3487 

1.0197 

0.7561 

3 

1.5855 

1.0256 

0.6469 

4 

1.7653 

1.0277 

0.5822 

5 

1.9105 

1.0285 

0.5383 

6 

2.0324 

1.0286 

0.5061 

7 

2.1375 

1.0284 

0.4811 

The  results  in  Table  E.3  resemble  those  in  Table  E.l.  In  both  cases  the 
value  of  E[zk]  increases  with  k,  whereas  the  standard  deviation  of  zk  increases 
more  slowly  with  k.  Thus,  the  coefficient  of  variation  decreases  with  k,  but 
not  so  rapidly  as  in  Table  E.2.  It  is  noted  that  when  the  mean  of  one  RV  in 
the  set  is  greater  than  the  constant  mean  of  all  others,  the  variance  of 
will  eventually  decrease  with  k  beyond  a  certain  point.  In  the  example  above, 
with  E[x.j]  =  0.79,  i  >  1,  the  maximum  variance  occurs  at  k  =  6. 

Computational  Errors 

The  results  in  Tables  E.l,  E.2,  and  E.3  are  numerically  exact  to  the  number 
of  significant  digits  displayed.  These  results  were  obtained  from  the  closed- 
form  solution  equations  using  double-precision  arithmetic. 

When  more  general  results  are  wanted,  it  is  convenient  to  use  a  numerical 


procedure  applicable  to  all  distributions  of  positive,  continuous  RV's.  The 
procedure  displayed  in  Annex  B  uses  equation  (12)  to  obtain  the  c.d.f.  of  z^. 
The  mean  and  variance  of  zk  are  obtained  using  (13),  (14),  and  (15).  Because 
the  procedure  involves  numerical  integration  at  each  recursive  step,  a  compu¬ 
tational  error  is  incurred.  With  rectangular  integration  an  integration  step 
size  of  0.002  is  judged  a  satisfactory  compromise  between  accuracy  and  speed. 
With  Simpson's  rule,  the  step  size  can  be  relaxed  to  0.005,  yielding  essential¬ 
ly  five  digit  accuracy  in  Etz^]  for  k  <  7.  Computational  errors  with  the  first 
order  procedure  (step  size  0.002)  are  displayed  in  Tables  E.4,  E.5,  and  E.6. 
These  results  are  regarded  as  representative  of  the  errors  to  be  encountered  in 
network  applications. 


TABLE  E.4 

ACCURACY  OF  A  NUMERICAL  METHOD  FOR  OBTAINING  THE 
MEAN  AND  STD  DEVIATION  OF  THE  MAXIMUM  OF  A  SET  OF  RV'S 

Case  A:  All  RV's  (x^)  are  exponentially 
distributed  with  rate  parameter  X..  =  1,  1  <_  i  <_  k. 


No.  RV's 


Exact  Values 


Numerical  Approx. 


TABLE  E.6 


TYPICAL  ERRORS  OF  THE  NUMERICAL  PROCEDURE  FOR  CALCULATING 
MEANS  AND  STANDARD  DEVIATIONS  OF  THE  MAXIMUM  OF  A  SET 
OF  POSITIVE  RANDOM  VARIABLES 


No.  RV's 
k 

Errors  for  Case  A 

Errors 

for  Case  B 

E[zk] 

Std  Dev  z^ 

E[zk] 

Std  Dev  z^ 

2 

-0.0004 

0.0000 

-0.0000 

3 

-0.0001 

-0.0007 

0.0000 

4 

-0.0001 

-0.0011 

0.0000 

5 

-0.0001 

-0.0013 

0.0000 

6 

-0.0016 

0.0000 

-0.0001 

It  is  noted  that  the  computational  errors  of  the  mean  value  and  standard 
deviation  in  these  examples  have  a  consistent  sign,  indicating  underestimation. 
Since  the  numerical  procedure  used  here  is  first-order,  the  errors  in  Table  E.6 
are  proportional  to  integration  step  size. 

Exact  distribution  functions  of  z^,  (2  £  k  1  6),  are  shown  in  Figure  E.l 
for  the  case  in  which  x. ,  1  £  i  £  k,  are  exponential  with  unity  mean.  It  is 
noted  that  substantial  positive  skewness  persists  with  increasing  k  up  to 
k  =  6,  In  Figure  E.2  is  displayed  the  c.d.f.  of  z^  for  the  case  in  which 
Aj  =  1  and  A..  =  1.2A^_j,  2  £  i  £  k.  Because  each  additional  (ordered)  RV  has 
a  smaller  mean  value  than  its  predecessors,  the  upper  tail  of  the  distribution 
of  z^  is  approaching  an  asymptotic  form  as  k  grows  indefinitely.  From  Figure 
E.2  it  is  apparent  why  the  standard  deviation  of  k  decreases  with  increasing  k. 
This  example  applies  to  those  networks  having  parallel  activities  in  which  a 


few  activities  have  much  greater  mean  completion  times  than  the  other  activi¬ 
ties.  In  this  instance  the  mean  value  of  the  completion  time  for  this  portion 
of  the  network  is  insensitive  to  the  addition  of  activities  beyond  a  certain 
point. 

Parametric  Analyses 

In  a  fully  general  situation  each  of  the  random  variables  in  the  set 
may  have  a  unique  functional  form.  In  such  a  setting  few  numerical  generali¬ 
zations  about  the  distribution  of  max  in  set  (z^)  can  be  made.  However,  one 
can  be  more  restrictive  with  respect  to  assumptions  with  some  beneficial  con¬ 
sequences.  Assume  that  the  functional  form  of  all  the  x^  is  the  same,  for 
example,  gamma  or  Weibull.  Also  assume  that  the  mean  value  of  all  the  RV's 
is  the  same.  With  these  restrictions  some  interesting  parametric  analyses  can 
be  performed.  One  analysis  of  interest  is  the  effect  of  distribution  shape  on 
the  mean  value  of  z^.  In  the  case  of  the  gamma  c.d.f.,  shape  is  affected  by 
only  one  parameter  —  8  —  in  the  function: 
f  X 

F(x)  =  BX(Xt)6‘1e“Xtdt/r(6)  .  (18) 

Jo 

As  8  increases,  with  fixed  mean,  the  distribution  becomes  less  variable  and 
less  positively  skewed.  In  this  case  the  coefficient  of  variation  is  l//8~. 

One  would  expect  that  ELz^]  would  decrease  with  increasing  8,  since  less  prob¬ 
ability  density  is  associated  with  large  values  of  x...  As  seen  in  Table  E.7, 
this  expectation  is  correct.  Since  our  general  numerical  method  was  used  here, 
results  are  presented  to  only  four  significant  digits.  Integration  step  size 
was  chosen  to  yield  four  digits  accuracy.  This  was  checked  against  exact 
results  for  the  exponential  case  (8  =  1).  A  generalization  from  Table  E.7  may 
be  of  interest.  Over  the  range  of  gamma  shape  parameter  shown,  there  is  an 
approximate  geometric  decrease  in  E[ z^ ]  -  1  with  increasing  8.  The  rate  of 
decrease  is  observed  to  be  greater  for  larger  values  of  k.  Clearly,  in  the 
case  of  zero  variance  as  8  approaches  infinity,  Efz^]  approaches  unity  for  all 
k. 
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In  the  case  where  F{x)  is  Weibull,  the  single  parameter  0  affects  shape: 

F(x)  =  1  -  exp[-(Ax)6]  .  (19) 

Results  for  this  case  are  shown  in  Table  E.8.  With  increasing  3  the  coeffi¬ 
cient  of  variation  uniformly  decreases  toward  zero,  as  in  the  case  where  x  is 
gamma.  However,  with  increasing  8  in  the  Weibull  case,  the  coefficient  of 
skewness  decreases  thru  positive  values,  becoming  negative  at  8  =  3.6.  As  8 
continues  to  increase,  the  coefficient  of  skewness  asymptotically  approaches 
-1.14.  This  behavior  (among  others)  distinguishes  the  Weibull  from  the  gamma 
c.d.f.  Coefficients  of  variation  and  of  skewness  for  the  gamma  and  Weibull 
distributions  are  tabulated  versus  shape  parameter  in  Table  E.9.  For  large 
values  of  6,  the  Weibull  is  not  really  comparable  with  the  gamma  distribution 
for  reasons  given  above.  However,  for  values  of  8  in  which  the  Weibull  is 
positively  skewed  and  has  the  same  coefficient  of  variation  as  a  given  gamma 
distribution,  certain  results  are  quite  similar.  For  6=1  both  distributions 
are  exponential  and  yield  identical  values  of  E[zk].  Consider  the  non-trivial 
case  in  which  the  coefficient  of  variation  is  1/2:  6(gamma)  =  4  and  8(Weibull) 
2.102.  In  this  case  the  form  of  the  distribution  has  little  effect  upon 
E[zk].  Results  are  nearly  the  same  --  within  ±0.02  --  for  2  £  k  £  6. 

The  c.d.f.  of  z^  is  displayed  in  Figure  E.3  for  several  values  of  k,  given 
that  F-(x)  is  Weibull  with  shape  parameter  of  2.  For  this  particular  distri¬ 
bution  the  mean  of  xi  is  located  at  the  54%  percentile.  Altho  the  skewness  of 
this  choice  of  F. (x)  is  rather  small  (0.63),  the  Gk(z)  functions  are  distinctly 
non-gaussian.  (Note  that  the  plots  are  made  on  Normal  probability  paper.)  It 
is  noted  in  passing  that  the  assumption  that  zk  is  Normal  is  generally  grossly 
wrong . 


A- 1 9 


TABLE  E . 7 

PARAMETRIC  ANALYSIS  FOR  THE  MEAN  VALUE  OF  THE  MAXIMUM 
OF  A  SET  OF  GAMMA  RV'S  WITH  SHAPE*  AS  A  PARAMETER 


E[xi]  =  1,  1  1  i  <  k 


RV's 

E[z^]  for  Gamma  Shape  Parameter: 

k 

1 

2 

3 

4 

5 

2 

1.500 

1.375 

1.312 

1.273 

1.246 

3 

1.833 

1.606 

1.498 

1.433 

1.387 

4 

2.083 

1.774 

1.631 

1.544 

1.486 

5 

2.283 

1.904 

1.733 

1.630 

1.561 

6 

2.450 

2.012 

1.816 

1.700 

1.622 

7 

2.593 

2.102 

1.886 

1.759 

1.673 

Distribution 

Type 


CV  for  Shape  Parameter: 


1 

2 

3 

4 

5 

Gamma 

1.0000 

0.7071 

0.5774 

0.4472 

Weibul 1 

1.0000 

0.5227 

0.3634 

0.2291 

TABLE  E.9B 

COEFFICIENT  OF  SKEWNESS  VERSUS  SHAPE  PARAMETER 

FOR  GAMMA  AND  WEIBULL  RANDOM  VARIABLES 

Distribution 

Yj  for  Shape  Parameter: 

Type 

12  3  4 

5 

Truncation  Effects 

It  was  stated  that  the  intended  application  of  the  above  analysis  is  to 
networking  problems  where  activity  times  are,  ultimately,  bounded  from  above. 
Thus,  the  reader  may  be  tempted  to  challenge  the  applicability  of  distributions 
of  x  which  yield  positive  but  unbounded  values.  All  of  the  above  methods  hold 
for  bounded  distributions.  Only  the  applications  emphasis  so  far  has  been  on 
unlimited  distributions.  This  situation  can  be  redressed  by  examining  trun¬ 
cated  versions  of  the  distributions  previously  considered.  Suppose  that  the 
RV  x  is  defined  on  the  finite  domain:  0  1  x  1  x  .  If  the  non-truncated  c.d.f. 
of  x  is  denoted  by  F(8,A  jX),  the  truncated  form  is  given  by 

F'(MU)  =  F(8,A',x)/F(8,A',xu)  , 

for  0  £  x  ^  x  .  If  the  untruncated  mean  value  is  unity  and  one  wishes  the  same 
value  of  the  mean  of  the  truncated  RV,  the  rate  parameter.  A,  must  be  adjusted 
to  A1  so  that 

fX 

u  F'fB.A'.xJdx  =  1  . 

o 

This  adjustment  to  preserve  the  mean  facilitates  a  comparison  between  results 
for  truncated  and  non-truncated  c.d.f.'s.  In  Annex  B  we  display  the  listing  of 
the  computer  program  which  implements  truncation  for  the  2-parameter  Weibull 
family.  In  the  case  where  the  shape  parameter  B  =  1,  the  exponential  distribu¬ 
tion  is  realized.  The  mean  values  of  the  maximum  of  a  set  of  truncated  expon¬ 
ential  RV's  were  calculated  for  several  values  of  x  .  The  upper  truncation 
point,  xu,  was  chosen  to  yield  convenient  values  of  F(B,A',xu)  such  as  0.990 
and  0.999,  etc.  Results  are  shown  in  Table  E.10.  These  may  be  compared  with 
the  untruncated  results  in  Table  E.l.  It  is  clear  that  truncation  of  the  xi 
values  has  a  much  greater  effect  on  the  standard  deviation  of  zfe  than  on  the 
mean  of  z^.  Specifically,  when  the  upper  1%  of  the  distribution  is  truncated, 
the  standard  deviation  of  ij  is  reduced  to  78%  of  its  untruncated  value,  whereas 
the  mean  of  z?  is  reduced  to  96%  of  its.  As  might  be  expected,  very  little 
difference  exists  between  truncated  and  untruncated  results  when  only  0.1%  of 
the  distribution  is  truncated. 


TABLE  E. 10 


a 


EFFECT  OF  TRUNCATION  OF  THE  EXPONENTIAL  DISTRIBUTION 
OF  x.  ON  THE  MEAN  AND  STD  DEV  OF  zk  =  MAX(x1>  ....  xk), 
GIVEN  E[xi]  =  1,  1  <  i  <  k 


Upper  Trunc.  Point*: 


4.8298 


6.9559 


Probability  Distribution  of  the  kth  Largest  of  n 


The  primary  focus  in  this  annex  is  on  the  statistics  of  the  largest  posi¬ 
tive  RV  in  a  set  of  k  RV's,  each  of  which  has  a  unique  c.d.f.  When  all  of  the 
RV's  in  the  set  are  from  the  same  distribution,  the  problem  addressed  here  is 
just  that  of  the  largest  order  statistic*  Order  statistics  pertain  to  the  kth 
ordered  RV  in  an  identically  distributed  set  of  n.  An  extensive  literature 
exists  relative  to  order  statistics.  Gumbel  ( 1958) [8] ,  for  example,  uses  this 
theory  to  develop  the  statistics  of  extremes.  As  pointed  out  by  Guenther 
( 1 977 ) [9] ,  the  c.d.f.  of  the  kth  order  statistic  from  a  set  of  n  RV's  with 
common  distribution  function  F(x)  is  given  by 


Gk,n^  =  ^(z)^*0  “  k  +  ^ 

where  Ix(a,b)  is  the  beta  distribution  of  x  with  parameters  a  and  b. 
distribution  is  expressed  as 

X 

Ua.b)  = - - -  ua-1(l  -  u)b_1du  , 

x  B(a,b)  o 

with 

6U,b)  ■  r(«)  r(b)  . 
r(a  +  b) 


(20) 

The  beta 


(21a) 


(21b) 


Guenther  indicates  that  it  is  computationally  convenient  to  obtain  Ix(a,b)  from 
an  equivalent  form  of  Fisher's  F-distribution.  A  numerical  method  for  evaluat¬ 
ing  the  F-distribution  is  given  on  p.  944  of  Abramowitz  and  Stegun  (1966)[10], 


[8]  Gumbel,  E.J.  Statistics  of  Extremes,  Columbia  Univ.  Press,  New  York, 
c.  1958. 

[9]  Guenther,  W.C.  "An  Easy  Method  for  Obtaining  Percentage  Points  of  Order 
Statistics,"  Technometrics,  Vol.  19,  No.  3,  pp.  319  -  321,  August  1977. 

[10]  Abramowitz,  M.  and  Stegun,  I.  Handbook  of  Mathematical  Functions,  AMS  55, 
Nat.  Bureau  of  Standards,  August  1966. 


Denote  the  upper  tail  probability  of  the  F-distribution  with  argument  y  and 
with  integer  degrees  of  freedom  parameters  Vj  and  v2  by  Q(y,\>j,\^).  Then,  the 
following  relationship  can  be  used  relating  Ix(a,b)  to  Q(y,v1,v2): 

Ua.b)  =  Q(..aU~.~  ,x?  ,  2b, 2a)  .  (22) 

x  bx 

Annex  B  includes  a  listing  of  the  computer  program  TEST.K  which  implements 
equations  (20),  (21),  and  (22)  to  calculate  the  distribution  of  the  kth  order 
statistic  for  a  set  of  n  generally-distributed  continuous  random  variables. 

The  distribution  function  F(x)  is  calculated  in  a  user-supplied  function  FUN. CDF. 
The  form  of  this  routine  illustrated  in  Annex  B  calculates  a  (optionally)  trun¬ 
cated  Weibull  distribution.  The  mean  and  standard  deviation  of  the  order 
statistics  are  evaluated  using  the  method  shown  in  equations  (16)  and  (17)  and 
employing  Simpson's  rule.  Despite  the  ease  of  implementation,  little  applic¬ 
ability  of  order  statistics  is  foreseen  to  networking.  This  is  due  to  the  fact 
that  parallel  activities  seldom  have  the  same  c.d.f.  When  they  (approximately) 
do,  it  is  quite  important  to  specify  whether  passage  thru  the  network  requires 
the  completion  of  all  of  the  activities  or  of  a  subset.  To  illustrate  this 


point,  consider  Figure  E.4.  A  comparison  is  made  between  Gg  7(z)  and  G7  7(z) 
when  F(x)  is  a  standardized  exponential  distribution.  It  is  seen  that  a  remark 
able  difference  exists  between  mean  values  of  the  completion  times  —  in  a  net¬ 
work  context  --  of  two  situations:  (a)  6  activities  of  7  must  be  complete  for 
passage  versus  (b)  all  7  activities  must  be  complete  for  passage.  This  illu¬ 
strates  a  suprisingly  high  sensitivity  of  project  completion  time  to  network 


ANNEX  B 


Computer  Source  Programs 

Two  MAIN  computer  programs  are  presented  in  this  annex*  The  first  (MAXG) 
can  be  used  to  calculate  the  mean,  standard  deviation,  and  probability  dis¬ 
tribution  function  of  the  largest  of  k  positive,  continuous  random  variables 
xi»  *  5  i  <  k,  each  having  its  own  c.d.f.  The  c.d.f.  of  x.  is  calculated  by 
the  function  FUN. CDF  which  accepts  as  arguments  the  shape  and  rate  parameters 
and  the  value  of  x^.  As  shown,  FUN. CDF  produces  probability  values  for 
(optionally)  truncated  2-parameter  Weibull  distributions.  Comment  code  is  also 
provided  for  calculating  the  gamma  distribution.  All  supporting  routines  and 
functions  are  supplied  in  this  listing. 

The  second  MAIN  program  (TEST.K)  calculates  the  distribution  function,  mean, 
and  standard  deviation  for  the  order  statistics  of  a  general  c.d.f.  This 
general  c.d.f.  is  calculated  in  the  user-supplied  function  FUN. CDF.  As  shown 
FUN. CDF  provides  values  for  truncated  Weibull  distributions,  just  as  this  func¬ 
tion  does  for  MAXG.  The  fact  that  FUN. CDF  is  identical  for  these  programs  is, 
of  course,  not  necessary.  All  utility  routines  are  provided  for  TEST.K.  Com¬ 
ment  statements  in  these  routines  explain  their  purpose  and  define  input  and 
output  arguments. 

These  programs  are  written  in  SIMSCRIPT  2.5  for  the  PRIME  750  minicomputer. 
However,  the  code  does  not  employ  features  unique  to  this  computer.  Cross 
reference  lists  are  included  with  program  statements  to  identify  variable  type, 
to  tabulate  the  locations  of  each  variable  in  the  programs,  and  to  facilitate 
the  conversion  of  programs  to  another  language.  Both  driver  programs  are  inter¬ 
active.  Program  input  is  read  from  the  terminal  and  output  is  displayed  at  the 
terminal.  No  external  files  are  used.  Since  the  output  may  be  lengthy,  it  is 
recommended  that  a  COMO  file  be  established  to  display  or  print  it. 


For  convenience  and  without  loss  of  generality,  the  random  variables  are 
considered  scaled  in  dimension  so  that  the  largest  mean  of  the  x^  is  unity. 

Inputs  to  MAXG  and  TEST.K  are  provided  in  response  to  prompting  messages  sent 
to  the  terminal.  Thus,  MAXG  requires  as  input  the  max  number  of  random  variables 
in  the  set,  an  indication  of  whether  these  are  defined  on  a  finite  or  semi¬ 
infinite  domain  (If  finite,  the  upper  truncation  point  must  be  specified.),  the 
shape  parameters  of  the  distribution  of  the  x^,  and  the  ratio  of  the  mean  x^  to 
the  mean  of  x^.  At  the  user's  option  the  values  of  G^z)  are  printed  out  at 
intervals  sufficient  to  permit  a  visually  smooth  point-to-point  plot.  Whether 
or  not  Gk(z)  is  printed,  the  program  provides  the  mean  value  and  standard  devia¬ 
tion  of  Zy  Zy  etc.  up  to  the  maximum  set  size  specified. 

In  the  program  TEST.K,  required  inputs  are:  (a)  the  cumulative  probability 
associated  with  the  upper  truncation  point  --  if  truncation  of  F(x)  is  chosen  — 
(b)  the  number  (n)  of  identically-distributed  RV's  in  the  set,  (c)  the  index 
(k)  of  the  order  statistic  wanted,  (d)  the  shape  parameter  of  F(x),  and  (e)  the 
rate  parameter  of  F(x).  Output  from  this  program  is  the  distribution,  G.,  (z), 

K  jl) 

of  the  selected  order  statistic  as  well  as  its  mean  and  standard  deviation. 


B-2 


tr 

Uj  UiUi  UlUlUiU/ 

Lm 

JJ  JJ. 

iOJ 

C: 

a.  aiaD  CD  qd  Ui  cl 

o 

3DDDDKD 

a 

ooco o?c 

OQQOOhD 

»“ 

* 

u  v 

<» 

vf  r- 

1  ►- 
J« 
ou» 
a  a  u.  u 

M  Z/L-  ►- 

a: 

aa 

tr  i  i 

a: 

ft* 

t-*q.  UlUi 

«i 

«*x 

«>  > 

u 

*  Z>Z> 

UI-JCO 

-»  u  u 

u 

>**«  <  M 

i 

'UJ 

>  r*  *3 

L. 

u 

*  *- 

u 

u  u 

«_* 

U 

_ 

-J-J 

Oj 

•-*  UjU* 

a 

OJ  Oj 

< 

H  »  >> 

UJ 

«* 

*x< 

»— « 

O  O 

1 

(g  </)«?  n*Oi*? 

u. 

a. 

a  a. 

«J  UD^“ 

«« ■** 

>  -JU.U4?  ^ 

Ui 

UJ 

> 

>> 

<£) 

o  o 

a 

u.  cu>  u 

ocu  oa» 

1 

CL  CL  CD 

u.  «t  at  _j  i.  ^ 

0333300 

►  >  ->ii  j 

VJ 

JOOCCJJ 

«*  Luu.  u. 

CULddClOO 

*— -j  a: 

< /) 

«  U!  <t  <  't 

c 

<t  IX  L  Ui 

-»  <V)QSOC 

Ob 

«  (/></>  ac  «x 

U  «J*ll  <x  < 

o 

Li  <1  h- 

jl.  o i  o'»  L/i 
«/  «t  K  •  ►  'i  i 

ur 

>  U.  ♦—  *u.u.u 
*t  Ui  •  *  UOOQJ 
a-  dul.  j^ua.Oi 

*  Ci*'  *a.<  •  «i 

•  a  z>u_  a  j?Z4 

aooir  D3u 

Ui  ►  u  o 

-J_>  CL 

< 

(Xi^UiUJUiUJuJUJ* 

a. 

a 

a 

•  i*. 

«■*  ►- 

UI  1  U  li  I»_  u.  1*.  I*  * 

•  >  '  V 

lie  UI  UJ  ui  UJ  Uj  loJ 

a  i'U.i  UU  L  UJ 

iJU  _J 

QoaDDa  *- 
UUUU  U  ►“O 


t 

l*< 

O 

a 


eg  a. 


«i  o 
uitd 


0 

CO 


c»  • 
e-J'/ 


o 

;* 

tt 

C 

to 


a 

x 


-*  u 
X 


3 

II  CO 


>  II 

<o 


DUJ  to  J<> 


«-«  > 
w  QC 


X 

U 

X 


u. 

D 


o 

z 

3 

a 


3 

tr 


>  •— « 
-i 


ii 

tt 

u> 


a 

>» 


to 

►H 

o 


a: 

UJ 


*-  o 
X  Q  «J  O 


*—  a  uj 

LCOU  2 


V  L 
L2 

z  — 


*-  *>4  O  J 


•  II  O' 

-a 


r>u  *»*- 
1  OUJUJU 

or  -j>~> 


u 

C- 

•1L7  Of 

u 

Wl- 

u 

.J 

►-  ft. 

*■< 

O  o  < 

X 

3 

X 

..1 

wt  > 

OiXU  UJU. 

H- 

_JU 

>- 

A<  U 

-ill 

'J 

-J<0  U  or  ►- 

O 

«J  *-* 

o 

too  Uj*m  f 

D  -B  II  II  II 

X  uj  t/>U_  LL. 

*-  X  «  CJ  (J  n  QC 

►“  <  o  o  ou 

l  —  c  o  o 

>-iXU  4  «  _JtO 

<1  L  •  -J  _J  L  S' 

-JD-Xuclll.  •  •«  Lj 

Q.C-I-  t/l  o  to  _J 

•*■(02  CO  2HhK 

<  CLUJ  3  Ui  40  3  UJ  3  UJ  to  «\J 
i-O  DJa  jVKJa 

4(0.2  too. 
3*-~D 
JXXIU. 

wcoa»-o 


to 


1- 

2 

UJ 

1- 

<2 

CO 

IX 

UJ 

1 

X 

Ui 

ID 

X 

£‘ 

e> 

h- 

D 

2 

It 

►— 

<M 

z 

•“ 

o 

4 

► 

1 

-J 

4 

1 

u 

*-**- 

co 

It 

u 

Ki 

> 

>- 

> 

W> 

x 

t 

o 

« 

U 

2 

NJ 

•-< 

Juj 

coco 

D 

a: 

ID 

»  • 

a 

w 

ID 

>> 

♦- 

> 

| 

aoc 

X 

1 

X 

1 

4 

>-> 

e 

1 

tl 

w 

OJ>- 

L 

l> 

X  » 

L- 

NO. 

z 

Uj 

<< 

CL 

o\ 

w 

ID 

oo 

X 

40 

a. 

I 

w 

• 

p 

oca 

— J 

♦  CVJ 

UJ 

F 

00 

t-J  ♦ 

X 

*o 

U.U 

• 

•  Uj 

K 

1 

w 

•H  • 

1 

1- 

l 

1 

►H 

UjCvJ 

W^ 

O 

1 

CL  4 

l< 

BE 

> 

4 

XX 

CO 

♦ 

K  1 

«x 

2  2 

wx 

2 

« 

UJ 

D 

UJW 

o  w 

o 

« 

IX- 

<*»>UJ 

•  <i 

*■* 

* 

1 

«MML 

UJ  • 

K 

4 

I 

— 4  '-O  ♦ 

h-U 

D 

4 

>  K-« 

«•**•  >  fVi 

UJ^- 

CO 

4 

1 

M>U.O  4 

JU 

►•g 

4 

1  1 

D> 

~<l-2  * 

Q _ J 

QC 

4 

IS  1 

M- 

>J^3X 

XX 

►- 

la  1 

I  -«»-  ii  ii  och- 

OX 

<o 

15  1 

1  *-*w 

H  m<UHO 

DO 

M 

4 

C  —U.O  II  02  *CO 

HO 

a 

4 

i  i 

D>  »<40)3mII 

211 

4 

UJ  1 

hhjkcocuq: 

WQC 

>- 

4 

H  1 

XffiKk  UJ**H>3< 

LI  W 

*— 

4 

H*  1 

WHO 

•  to  J*  > 

LX  > 

* 

la  1 

X  i  coo  — J 

to 

-I 

M 

CO 

«1 

tr 

o 


<i  g~i  ►— 


a  <r 1 

41 

t  >  • 


^CgK;  «y 

[■  si> 


u  o*o^*aoff'’  <t if  vflr^ccer x 

^  KN h» SNN  ^  NS  4 


> 

O 


ix 

la 

|w 

IQ. 

I 

UJ 

IO. 

Iw 

IX 

H- 

I 


I  4H2HHHJKH 
IO  UJ  O  D  UJ  UJ  UJ  «  UJ  3  Li  LI  CO 

. II 


rx.  co<r  -j_j_JD_jac 
I-  Uj 

x 

in  i 


III  LULL. 
^-4_J-J»~' 


> 

I- 

OJ 


♦  ^  J 

«  QC 

2  W  UJ 

«S  >  *-*2 

UJ4JZ  II 

IU.<UJ 

I  «u< 
0C»“*  Jr< 
<crsu 
><30  •  ►- 

II  tO COD  ZOO 
QC  l»  II  2*-*D 

COWO>D<KX 
j>->woaai- 
<  >- 
*i-*-i- 

-JLJ  UJUJU. 


g 

12 


K  _tc\iK  <r  ir  v^r-  cr<T 
locacccocfrxooooaa. 
I 
i 


C  w(  j»n 
O'ffff'O' 


<j  iT  \T  a  ^o»h 

(J\  (MJs  O' f  t  4 


B-5 


■r  c«c’  yrt'crPt^T  11.5  p-r  ppi**e  syste“s,  pelepse  ?. 

pPi^ESCE.  TD.S  LoOK  .xREr  »KOf  YPLTST  «  T<  s CE  ’  !r  CEC  190P 


1  ►- 

> 

> 

1  </> 

0 

■ 

• 

1  a 

c 

0 

1 

• 

c 

CO 

1  «* 

•-* 

2 

• 

1  X 

<1 

i 

I  X 

II 

a: 

<1 

L_ 

1  0 

k 

UJ 

a 

1 

> 

• 

1  a 

a 

*—> 

1  O 

0 

t- 

*- 

1  u 

u 

»— 

1 

00 

> 

1  <Vi 

«•* 

c 

*a 

1  * 

■ 

a 

* * 

1  • 

• 

c. 

<~ 

*• 

•-< 

lb 

*—■ 

u. 

1 

K 

w 

1  > 

0 

2 

• 

> 

1  r*. 

»► 

** 

u 

H- 

1  u. 

• 

IS* 

> 

c 

0 

Oi 

1  « 

u 

0 

u. 

• 

1  ♦ 

• 

UJ 

•X  0 

r*M 

X 

«— < 

1  Cm 

0 

0 

fM  O 

Hi 

t  * 

• 

• 

•  X 

u. 

X 

X 

1  « 

0 

<\l 

CM  I 

O 

»-« 

H- 

a. 

a.  ~ 

X 

1 

0 

• 

•  • 

■ 

«4 

<• 

1 

2 

H 

u 

X 

1  > 

4 

a 

0.  > 

• 

UJ 

IP 

1  < 

VU.U. 

a 

Ul 

•  • 

1  -J 

• 

u  *»• 

u.  OQ 

• 

X 

K- 

*  :»</> 

I  V 

U 

a-c 

aou 

<_> 

»- 

-i 

•  a:  3 

j  a 

• 

o>o 

0  11  x 

CO 

*  X 

j  H 

O 

•u. 

fsj  »«  II 

Of 

UJ 

u 

rl 

}  Ow 

• 

2C2J?2p< 

0 

2 

0 

•  -1 

|  9>U 

Ci 

oaxwD*? 

u. 

*-» 

<*J  H- 

•  <14 

|  *-«b-o 

U  i»  <OU.>w? 

w 

«/)*'<// 

L*J 

n  a-o 

L 

I'  *»  > 

*  11  u.  > 

(/ 

(/  - 

0 

—  ►O 

CL  2 

I  ~u  • 

X 

2  U.OU.K 

O 

U.  II  <X) 

■  Of  X 

<PuJ* 

I  is* 

*- 

•n^oiiououj 

2 

ZU- 

l~ 

oc  a  ►- 

>->-! 

|  ^Hli  j 

^  >K  fsiO  *  Q.  > 

hOH* 

3 

u 

I  >  NlUJ 

UJ 

>14. 

O 

a 

-JU2COOJ 

X 

a*  u 

n>jo 

fs 

U  OC-w 

U- 

«* 

►- 

-j» 

i  UMUill 

M 

OOII  UUUUi* 

UJ 

cm«(qc  Ik 

0 

«  «x  ►- 

1  OON 

-4 

1X2JJJJ 

X 

jahu 

*  Cl  - 

1 

« 

a. 

au- 

CO 

*  0  ►- 

I  t-t-t-t- 

♦-►-a: 

0 

►- 

M 

•  0  CL 

1  UJLJUIUJ 

►- 

U<U-0 

0 

3. 

XU 

c. 

•  ja 

1 

»-« 

-J_|U 

—t 

►— « 

00 

• 

1 

a: 

>■ 

1 

► 

a 

*— 

i  *  * 

•  • 

• 

•  • 

•-< 

rg#o<r  it  si  r 


I  CD 

'  «j  r-  1  ooJT' ♦  ip  v£f»a  <r  0^04*0 «-m  \or*oro  « 


I  C.  IUV  < 

|  O  ~  II  II 

|  ZCNUi.  u. 

|  ~  2  *  *U  UJ 

|  C  *N2C  c 

1  *-*«»-*  r.«x  Jvouu 

|  (\.C\i>  0 C  •  2  .JU  4/) 

|  •  •  «C  O 

|  >>HC  rHhoilQUJJUM 

I  HOH  II  IIIHKO«J*J>“ 
I  0^111X13  2  u-  < 

I  II  il  tlDDDWHh  x  a 

fr/> «-* CM *-C/JC/> COX  II  UJU  H-  _j 

b-Q.0.3OXX2«Ji-i  C  < 


JU/WU-uJUJUJUjO 
ijw  Jjj^^U. 


1 1  lo^CMm^ir^cr-flDir^-^cMiojr  *n<i)C  «.cr 
D^*  tototomtototoromto*-  +  +  <t  +■  ■+  +  +  4- 


U  U. 

*-  **Q 

x  isjo 

2  *  » 

a  i^o 

u  a  « 


2  u.  2 
*-•  Q 
«x  u> 
*-  «u. 

CD  20 

o  no 
u  x 
O  II  II 

I-  *  CL 

U  _J 
J  QO 

<  OX 
uc 

O  I-  h- 
UJ  UiUi 


'■Dm 


C  I  N 

a.  x  a 

•u.  at  x 

i'O  a 

DO  vOiT  a 
u  «  *<a  ji 

«  D  Mt 

^oo  a  uj 

-  —  >  X  O  • 

>u.  a  •  a 

>liOO  _J 

a  auHacoH 

QU>  UM  ll«- 

o  >  i  u.  a.  au. 

X  ti  OULOOH 


D  >  <  «  > 

V)  *  I  JC>N 
<Ox»  •  «UJ 
♦  •CNJKIKX 

•H4 


U  U- 

v>  | 

uu 

iu  | 

oo 

• 

-j  1 

X  X 

CVJ 

CD  | 

«*  I 

XI 

A 

**  ! 

a-  1 

aau_a 

X 

«  1 

JJOO 

»  1 

OL 

I 

1 

i  Z 

a 

X  1 

a 

o  oh-  a 

M 

O  1 

X 

*-  HUU 

O  1 

D 

aa«« 

2  1 

00 

«X4XX 

• 

«  1 

X 

a 

ac  1 

X 

oo 

X 

w 

LJUiU-U. 

w 

Ui  1 

L> 

a  0.00 

Kj 

>  I 

• 

XX 

a.  I 

a 

LLZZ 

• 

*-  l 

cr 

OO 

V) 

a  1 

cs 

UjUJaa 

Ui 

00 

co 

>>aa 

-i 

c  i 

c  * 

»"»►  «« 

o 

• 

a  I 

*  X 

VjujOCQC 

«A 

a 

I 

♦  D 

ac  aou 

a 

X 

«  1 

X 

=>3aa 

a 

** 

«  1 

DX 

UUHK 

ki 

</>  * 

aazz 

> 

a  | 

XCJ 

<*>0Caa 

a 

o  1 

1  VI 

X 

o 

X  < 

>■  >->“>- 

c 

X  I 

=>  • 

CD  OD  0)01 

o 

• 

3  1 

v>~ 

u. 

X  t 

xc\, 

oooo 

•a 

• 

H  | 

X  — 

UJULJU 

a: 

a 

X  X  | 

JJJJ 
DDD  D 
UUOU 
JJJJ 


X 

JX 

II  a  •Q.3WUZW 

•  o  i 

1 

>H>UU|H 

OOOO 

< 

ac 

a 

a 

l 

o 

o  • 

o 

-/•Hfi.  «  *~D 

•  -»-i  aa 

1 

N^^OOJ 

X 

e 

2 

o 

i  ~ 

tf  1 

X  f'4 

mZ 

—  u  Drsii:  «c  uc  X 

•  v  yr  X 

1 

u>a  ♦  * 

*  *  •  ♦ 

a 

a 

DOO 

»  iC 

*  • 

w>a:  «  ♦  -Jli- 

•  >-  ^  O  Ui  >•’  a 

1 

«*rs.  •>  t  i/i 

«  •  «  « 

a 

o 

-J 

2 

DZ 

1  w 

>  u  oa  a  a  u. 

•  <  >>-</)  Ui 

1 

II  >^3  »u 

«  *  *  ♦ 

X 

o 

U) 

(A  - 

c 

XUJ 

1 

• 

i.JU> 

• 

a  aQ.  jj  • 

•  -»  a  ^-ctar 

1 

H 

«  «  «  « 

a 

o 

a  h  cp 

a 

IQ 

I  ^-vQ. 

D 

f  o 

a 

ooaooau. 

j^QoayH 

1 

N'OX  II*" 

•  •  •  • 

X 

S’ 

2a 

I  >>a 

a 

oo 

• 

a  >  dooo*- 

<>*UO«>J 

1 

>-MH  III  -1 

*  *  *  • 

QC 

< 

a 

-<oa 

3 

1  a« 

• 

o 

X 

4  x  ao 

1 

>«-  II  X.D 

*  «  •  #vo 

O 

X 

a 

JUZl/lffi 

1  IQJH 

oo 

• 

a  a  aQQ 

31  X»  <Vl 

1 

*i»OD(rt<r 

a 

a 

«* 

Uh3h 

1  ii  M  a 

o 

CO 

a 

uiUiLjooa 

*  JU. 

1 

a  Kuo  oo  x 

* 

O 

a 

KXCCX  cc 

icoacja 

s 

■a.  "o 

*  -*  *-*  o  ►- 

! 

«i><xxa 

lllh 

(O 

a 

X 

aaaa 

ms 

u. 

*  LZ 

i 

II  P 

a 

a  a 

CO 

K 

J 

a 

« 

Mil 

a 

K) 

a 

a 

a 

oa 

t»aaa 

3 

«3 

*  oo 

i 

yuuuuoa 

04 

D 

II 

X 

ocu 

o 

DC/5 

Maaw 

Uu> 

o 

•  oo 

1 

J  JJ_j  jO*  h- 

AAA 

<o 

X 

a 

00  a 

JC« 

KJJJ 

« 

«  » 

1 

wojcucvja 

a 

ac 

>- 

a  o 

1 

»• » -« cj *o  ir>  oa  at  cr- *<M»o  *  ♦  u"> aiafloiT'O 
w>irnra,!u  u'u'>W'«r  tr  xi sc  >0^;  *  ^ v£ 'Dn/i  wj'Xk 


aoi»o  #■  mtfr"* 
aaaaaKK 


•  »  •  i 

(/)•  »  »  I 

U->lL> 

ouoy 

Q  O®0*t,xi 

Z  Z  *•  k-cl-cd 
X  Vlttfl 


<^W«ifHCN®0'4’ 
at  cc  ao  co  <30  X  ao  00  dD 
aa*-aaaaaO 
ac 
& 


B-7 


?F-T*cvs  -  rrCl,fNCt,ICt'U'  CHK.xOf 


frGlTNCE*lC*$UhCbK*XREp  »  NC  E  L  !  $  T  ,  Tc,rr’  '  "  !*  DEC  19M»  1  :4o;u 


* 

If  u 

♦ 

♦ 

O' 

*-* 

CM  if' 

♦ 

sA 

♦  « 

moo 

« 

IT 

K)  CM 

m 

am 

cm  r- 

CL 

o 

o 

*A^ 

3 

O 

♦  ®  o 

CM 

mm 

m 

O^o 

CM 

CM 

CM 

3<M 

CM 

CM 

rl 

CM 

H 

•<H 

u 

o 


u 


m>  m 

mo 

CM3 

vO 

a 

mo' 

h- 

H 

m 

trm 

sfl 

oco 

f- 

CMft 

mo' 

o 

m 

mo 

« 

m 

♦  sACb 

O' 

♦  o 

m 

me 

H 

CM 

*4 

CM 

3 

HM 

3 

3 

*4 

H 

r<>  i 

« 

«  ♦  ♦ 

• 

* 

K> 

r- 

r-r^ 

m 

®m 

O'O'  « 

■-4 

O 

CM  O'*-* 

CM 

NffiO 

♦ 

r-® 

r* 

M 

©♦ 

♦m 

o 

if*3 

mm  r- 

sC 

m 

♦  mo 

O' 

a  ♦a 

CM 

HIT 

CM 

H 

CM 

3  CM 

•H  *H 

*“* 

*H 

H 

3 

</j 

UJ 

* 

♦ 

♦ 

o 

—4 

sA 

•>ft*A 

CMm 

sA 

3(M 

®mo>  b 

r- 

If’ 

m 

3 

mm*-* 

m 

O'  sA 

♦ 

2 

Sfi 

(OITC 

♦  *-' 

m 

in*-* 

HNCJIO 

CM 

U 

♦ 

♦  mo 

O' 

00  ♦O' 

CM 

HffN 

*H 

UJ 

3 

CM*-* 

•— i 

H 

•HCM 

hh(V 

H 

CM 

3 

3 

CM 

3 


•ti 

Cl  a  A-  3-«ft  OO 

«  cm  cM3ft  ♦  ‘ 
3  r-i  CM3  3 

C 


in  oac^ioui^actf'  c 
US  me  O'  C>*-»f*-CuCM  3 
r-*  «HCM  *HrH»-«CM 


aj  o  b-iO'A 
♦  ♦  mm® 


out 

®m 

•-•CM 


3  CM  CD 

♦  ♦O' 


♦ 

CM  <JL*  O 

CM  r-»(T  O 

CM 


CO 

a: 

o 

CD 

X 

ID 


*  m  + 

co  ojoo  c  hocnj 

<Mro<M  f~  iT<r  CMibojfO 
CM  33  3 


CT>  K)h-OO*H(MVDf-CC00 

♦  COflNC  3f»CM3^  s£ 

3  «HCM  3**3<M 


O'  ctsacc  vLirm 
h'  ♦  cd  h  j  ♦  ao 
*hCM  CM 


CM!"**  O  CM  Si.  -c  in 
CCM  C^CM  ♦CM® 
CM  HCNJ 


HCMOvOCO  O 
CM  ft~>  ♦  3  ®  sA 


£  mr>  foco*o^a'f^mcDa'»-'u>®o*-«#ocM®o«H*niOh>-  ♦  cm  mme\if-  ♦«♦*>♦♦ 
*-CM*-.r  MOf't  — *a«-«m  •-^♦•-<  ♦  C.  ♦  ON  _<  f*.  cm  3  ♦  *►*  CM  nmiOCH^ODO  «-••-* 
j  Cs.  rH  «hCM*-h<m  HH(M  CM33  CM  3 


•a  ♦  r*-  O'  m  m  Cb  <m  cm  ♦  ♦  o  cm  ♦  r*  m  © 

CM  3  vA  3  ♦  3|(  Cb  m  ~  3  CD  CM  3®  CM  ♦ 
CM  3<M  3  « 


(X 

a. 

QL  QL 

ft 

ft  ft 

ftft  ftft 

UJ 

UJUJUJ 

UJ 

UJ  UJUJUJ 

UJ 

UJUJUJ 

33 

333 

3 

33 

3 

33 

333 

33333 

UJ 

-J 

•J-im. J 

-J 

JJJJ 

-J 

3  33 

30 

333 

O 

OO 

O 

33 

300 

00300 

Cl 

a. 

►— 

(DffiCD 

m 

Q.CL  CLiCD 

CD 

rnmes 

*»0Q3 

CD  CO  CD 

3 

33 

3 

a)  cd 

(DU  3 

330333 

o 

3 

X 

3  3D 

3 

3333 

3 

333 

X3»- 

333 

►- 

b»- 

3 

33 

3H- 

33^33 

I 

C 

UJ 

OOO 

O 

oooo 

C 

OOO 

LJOZ 

OOO 

Z 

2.  Z 

Z 

OO 

022 

ZZ02.Z 

O 

>~- 

OOD 

o 

OCCQ 

O 

ooa 

KO*-* 

ooa 

3 

33 

3 

OO 

033 

33033 

a 

o 

1 

o 

1 

a 

i 

3 

w 

w 

w 

W 

u 

c  >*Am 

CM 

K  sj  mMJ 

♦ 

CM  0 

*H 

♦ 

mu 

b 

u  -  a 

a  o- 

C-3 

a>cm 

•— • 

sA  HN 

mm 

m 

m 

m3 

3(f) 

o 

O 

OOO 

Q 

O  CjO 

o 

oo 

DDu 

o 

o 

oo 

o 

oo 

OO 

OO 

a 

a 

cr  era: 

QL 

era:  oca: 

QL 

OLCC 

QLCLCC 

a 

ex 

CLOC 

a: 

ex® 

ft  ft 

ftft 

O 

o 

otc 

C 

C  QL  cc 

o 

CO 

ccc 

o 

o 

CC 

o 

oo 

cc 

cc 

U 

* 

u 

jaa 

3 

303 

3 

33 

333 

3 

3 

33 

3 

33 

33 

S3 

o 

UJ 

Uj 

U<  UUJ 

UJ 

U>  uJUJ 

UJ 

U  UJ 

33  3 

3 

3 

33 

3 

U>3 

33 

33 

-J 

-1 

O-J^J 

mJ 

— J  Jj 

3 

33 

JJu 

3 

3 

33 

3 

33 

33 

33 

Ui 

c 

a 

a.  a- a. 

ai 

a  c  a. 

a> 

a  gd 

CCC  a 

CD 

er 

a  CL 

QL 

cd  a  ‘ 

Oft 

ftft 

•u 

HI 

<1  «*  *Q 

«» 

♦  u  «« 

HI 

HI  «c 

♦  ♦  ♦ 

HI 

HI  « 

«i 

ftft 

ftft 

u. 

*- 

►  • 

►  ►  .►H 

•-’ 

•-«#-« 

H 

►« 

►  3 

33 

33 

a 

a 

a  aa 

QL. 

U  Uj  a  (X 

a 

a  a. 

era  a 

a 

a 

ft  a 

a 

ft® 

ftft 

ftft 

-J 

♦ 

*a 

H»  ♦  ♦ 

♦ 

<<<< 

♦ 

♦  ♦ 

^H*  4 

HC 

♦ 

<  ♦ 

♦ 

ft  ft 

ft  ft 

ftft 

> 

> 

>>> 

> 

>H» 

> 

>> 

>>  > 

» 

> 

>> 

> 

>> 

>> 

>> 

u 

u 

U- 

UUb 

Uj 

UJ  <  oo 

UJ 

tJUJ 

333 

UJ 

u 

3  3 

u. 

33 

33 

33 

UJ 

UJ 

> 

> 

>>> 

> 

>^» 

> 

>  > 

>>  > 

> 

> 

>> 

> 

♦  > 

>> 

>> 

ft 

3 

*-*►-«*— 

M 

*— •  i  i»  i 

►H 

a  »-••-* 

•— >-*  *-* 

33m 

M 

»<!■  i 

3 

33 

UHm 

33333 

a 

> 

to 

V) 

<nv7oo 

O) 

o>  jmm 

00 

JJOOCO 

(/)V)V) 

22(0 

00 

m«/> 

00 

(/>(/> 

zoom 

zmeozz 

¥- 

at 

U 

cl  aa 

a: 

a -a  era. 

a: 

t-naca: 

a.  oca: 

•hohQC 

cc 

QCQC 

ac  33  ft  ft 

3a:  * 

3ftft«^3 

i 

3 

D 

333 

3 

3®  33 

3 

K33 

333 

»-i-3 

3 

33 

33333 

K33 

33033 

o 

UJ 

UUb 

U 

ocoo 

o 

3UU 

U3U 

330 

O 

oo 

OflDOOO 

300 

30033 

00 

3 

'*j 

uJUJUi 

UJ 

UJ-JUJUJ 

UJ 

033 

33  UJ 

003 

3 

33 

3ftft  33 

033 

C  U  300 

Cl 

aaac 

a 

O.OQ.CC 

CL 

a  a  a 

a  a  a 

a  a  a 

CL 

ft  ft 

ft  3  3  ft  ft 

ftftft 

ft  ft  ftftft 

(/> 


C 

a. 


•  3 

ui  a 


1/ 

3 

3  U. 

33 

3 

u 

a 

3 

3  a 

OO 

ft 

• 

o 

a 

U- 

a 

3 

•  u 

Oft 

O 

3 

u  uu 

*- 

*4 

J 

a  < 

•  >  X 

3 

ft  Nl 

3  O 

•  «o 

CD 

• 

•zz  •  • 

»— 

•Si 

i/ 

Ci 

> 

u.  u.  u  U. 

3 

X  3 

00  X  ft 

ZZ3 

a  > 

o 

X 

x  Xft  ZO 

J 

3 

z 

tO  403 

ooa  a 

O 

n>3 

33  3 

330 

*  CM  ft  ft 

o 

ft 

«  ft  330 

*1 

*4 

•J  ft  CL. 

a. 

oooo 

O 

OOO 

03  U 

U  3  X  3 

3ft  - 

J 3 333 

JB  X 

xxa  ax 

d  n 


c 


UJ  a  «co 

S3  H  H  Hrtrt 


u  c 

• 

•-1 

* 

• 

V/) 

u.  CO 

» 

«. J  ►-» 

L 

a 

►* 

CO 

•-•-tr 

X 

a  uj 

</ 

<1 

>3 

u. 

L 

7 

J  a 

o  «* 

a 

c 

u. 

-  > 

«4  U 

tr 

a 

c 

u 

u 

U<3 

>  t- 

• 

U-. 

►- 

-1 

0 

» 

V- 

«v 

* 

r-> 

>-  • 

CL 

u 

► 

1/  c 

O 

a 

oa 

-J 

o 

X* 

a 

o 

a 

U> 

•— * 

♦£ 

UJ 

u. 

QC 

•-H 

U> 

u 

«/  Hv 

CL*- 

* 

n> 

u.a 

CL 

♦ 

*-'C 

* 

h- 

«/■ «- 

luU 

-> 

c 

u 

c  o 

U- 

*->X 

<1 

UJ 

<1 

a  • 

<1 

«_• 

a 

2 

< 

I  ^ 

u  1/ 

u  o 

U  2 

w  a. 

c  *--> 

• 

u  * 

•J 

K- 

CL 

JJ  ♦ 

•  * 

a 

•  3 

• 

■«  O 

a  -* 

* 

X  .x 

U- 

a 

X  o 

t.j 

•  *-« 

a 

u.  —  «j 

uj  a 

o 

aa 

O 

<  ^ 

1  < 

• 

•  »— 

O  w 

► 

K 

O  l/l 

• 

a  u- 

•  1 

►>» 

o 

U  W  (/■ 

•-•  w 

UJ 

• 

UJ  o 

< 

UJ  X  ►-» 

w  U. 

a 

kM 

X 

-J 

0*0 

• 

X 

ti 

►-U 

U- 

UJ  LJ  • 

X  UJ  CL 

• 

I 

• 

*- 1-  a 

X-  •  X 

3T 

• 

cov- 

o 

73  7 

••  *-4  UJ 

X 

—t 

UJ 

2 

I-.OC  -C 

3  1 

o 

a 

♦—  u 

3 

H-  _J 

X  o 

*» 

«1  l 

u 

i  •  a  ii  a 

2  X  • 

a 

M> 

<X  *“  UJ  C  ►-  JC  LkJ  •“* 

Vk 

_  CO 

'-v~  «1 3 

• 

o 

oa 

UJ 

t/)U-  *  X  -J 

»-*-•  2 

i_i 

-JuJ 

UJ 

*4  »ftiHU  VU  k~ 

• 

<  K- 

-1 

XCl  *-•  • 

a  a  uc  uj*"U. 

• 

UUi 

aj 

-J  **  II  Ji<JH 

LOW  JO 

UJ 

> 

< 

X.O  7 

O 

a 

u 

?  «* 

>—« 

UJ  II  (L23U. 

hLjJZ  • 

«* 

c:  a 

ac 

22^cxac»- 

oactx  a  7 

tl 

i. 

►-  «i 

u 

DH 

W33 

3 

<  / 

H-a 

> 

XI- u. 

L3 

► 

o 

UJ  U<uJ  UiU 

X-UJ» 

U 

K- 

7  U1 

Uj 

OJJKM 

o  <r  • 

i/» 

V 

aa 

1 

. 

«  H- 

i 

c 

•  •  • 

*  • 

•  • 

U-  * 

*  »  • 

*  » 

»  » 

Uj 

(// 

c. 

—  c 

x/  f-  a  <r  t  «-t  cm  k>  ♦  in  «a>  x  cu<r 

*-  ' 

*■’  ^  H  •— <  H 

2C  O  ♦  ♦  ^C4fO 


UJ  1. 1 1. 1  Ui  1.11.  Il.li.ll.fi.  I 

O  CD0D(DC&UJCOaQCDUJ<X> 

o  3333h333h3 
a.  0000200010 

OOOOhOCjOHQ 


o-« 

C\J  K-.  v£  OJ 

a 

O 

oca: 

a  •  •  *a;ac 

Oft 

COCOftft 

Ui 

31 X 

32ZZ« 

u 

7 

UJ 

UJ 

3 

-J 

L 

00 

a> 

«Uj 

«  U.Uj 

a. 

^  -« -J 

a  a 

a  cx  co 

UJ 

x  < 

< 

»►* 

> 

U. 

a; 

OCtL 

U-< 

Ui  <<t 

UJ 

Uj  >> 

ql  *-  ujUi^zzae 

X  W  JZ2  C/)  UJ  UJ  Uj  _J  _J 

►-  oc«H^ai:z 

303I“>“3  33  30DCD 
0033000000 
UJOOy^attjj 
ttoa:ixa<  i  <oo 


u>  u. 

Li  C3  • 

*  au.u  u 

4  •  •  •  2 

cu,a?  so. 

KQK3  ^  CM  Kh 

«ou'u  s  aa 


B-ll 


FUVCTIC*  FU‘.FOF  <  F 1 »  T,  DT> 

1 1 

••'UNCTION  CALCULATES  THE  DIFFERENTIAL  ro-'PfPILm  IN  TnF  INTEF VAt  t  <T-dT,T> 


</ 


LJ 

OL 

UJ 


uJ 

a 

u. 

o 

CO 
OC 
U i 
CD 
X 
3 
z 


a 


u 


*  * 

00  X  IT> 


HCOHHMH 


UiUJldLjLJUl 
Ui  -J-i.J_J.J-f 
c.  CD  CD  CD  CD  CD  CD 
O  33D333 
2  OOOOOO 
DDQCQO 


U 


eg 

u 

• 

*— • 

u 


u. 

c 

c 


Ui 

o 


a. 


•  •  •  • 

O  c  OC 

z  zzz 


u 


u 


ZU UZ2Z 
uiZ  ZUJLiUJ 

S  HHtU 

033000 

a-coaaa. 


>  u  o 

•  a 

o  u  <r 

• 

T  Z  O 

F  3 

*“  u 

• 

X 

Z  H 

>-  O-UIF-U. 

DJt/i  -*Q 
h-  *-«  Q- 
VU3Z  • 
a. a  xZ 
►-  u  3  3 
i  h-  u 
U  ►“ 

«  o  a.® 


u- 

•  z 

»  Ui 


uu 

ao 

ud 

•  • 

2  r 

»-  33*H(N, 


rK><rrfv/f--<r 


F  C(iUENCC  CHK  ,x  9EF  .M59  *Dl  T  TT  .  T<-ACE  3 


r- 

CM 


4K>  X 

K>  4  4-  #»  > 


l  4\J  D  K)  CDOUl 
»OC'4  4CM«Vtf  »0<\ICM 


4  •‘•cg»nr-oo0'4r><rf,,“tfi«- 

f<  /  f  ' •— <  4  »— «  *-•  CV  HHr 


UJUJU  UJU) 

UJ  O - IJJ 

C.  UJ  C  OJOjOj 
O 

a  « ooc  o 
«oooa 


*  i 

CM4"  h- 
CM4  4 


CM  4  IP  4 


C0< 

rt^^K) 


4K><Q4 
«■**■>  Ki 


it  OhKiiohuhK) 
K;^K)KjM»-'fO  K3 


K« 

UJUJUJ  UJ  UJ 
003  «J  3 
UJU,  03  £D  QQ 

♦“*-!>  D3 
220  OO 
►*i-Q  OO 


a  a 
ct  •<* 
o  oo 

3i  Z3 


</  »■ 
a  »• 


OL  UJ  UJ 

u.  >>>  > 

a.  u»L-i^*^*-*  ujUj»— 

>-  i  i  2i-  «/> 

•-«*-  .ucac  u_i333  3 -j  3  3  3  3 

*-»-  3  3  3UJUjUJUJUJUJUJUJUJUJ»-*-3 
3  3  UU  CJtfi  CD  CO  03  00 CD  CL  OO  QQ  3  3  O 
i  <OUJUJUJ«J  ««•  «J<«3«a  1«J*CCUJ 
auaaa  jjjjjj^jjjaira 


U  3  -f  •  O 

•  a  x  i,(3o  J7 

*  3  r  a  '  WCM^IT’-  r  r  -f  3 
nr.  o  u  <  ^  m-4^-4 ~nf  ^  oOuj  or 

*’  Uu  OO  JJ-JJ-JJ-J 33  JH.K  x 


B-15 


»-*  w  I 

O  L  u  *' 


</5  o  -J  wu 

Dr  w  or 

x  d  »  li 

cr  w  ►-  qc 

►-  u.  ►- 

Lj  •  •  -J  • 

l-HUi  >  LJ 

•  -J  3  u  -J  -« 


3  ft  f 
-  0«T 

>  C  COuj 

'  II  X _J 
*•  Cfc  I. 

I  J  «x  a_  x 
5  II  OJ’-'fft 
v  x_ia 


r  ►-  lx  n  -jr  *-  u. 

HDf'MU.  C_ 

^OkiOU-itOK 

a.  xv)crUQQ  cr  c/> 

^  o 

Oil-  ^ 

I-  KUiU.  ►“ 


XU  I-  * 

4  •(/)  * 

I0.HO3  a 
X  XI—  •  *  .1 

♦  U-  --  «-*Q  Li 

4  *  ♦  >v  ► 

w  r  </  <i  t .  i 

Oj*— <  *  f—  •  CL  *■ 

a.  -J*  u-^<  ii  a 

u.<  ii  oi  I.  a  « 

WJIC'-M  O  - 
a  oti  a 
a  ii  jo  ii  lju.'  - 
or  DUJi“  0.1-0  0 

UJ 

Ihl-HKhl-  ► 
KLiU  M  LUL  L 

U  J  JL  _ i _ ( _ I  _ 


r 

ru  v 

♦  c.  ~ 
uvr  a. 
w  ^  tr 

•  k.>u 

OHD 

r  ^  to  o 
ac  u«  2*-« 

lcc  «  o 4 -i 

u  C  II  JH 
*-  a  ii  .njcj 

*'  .“c  j  a  ► 
«*»  -o:a 
x  c_  juju  a  c 

IX 

>  a 

X  «“►-  >-o 
D  nut ou 
ooa  _i«j< *i 


o  *■*  w 

ID> 
1  S. 
D 

a  o 

u  a.  '■»  •  ♦ 

I/1  Ll- 

-0.0X0 

awuu< 

0->>llt. 
UOCjH  II 
X*  •  <_>  X 


u  a 

oo»-i- 
u  out 
— I-I-J 


<t  a  \inu!j  •  a  c  in  o.r-^ 

0-  Q. 


4 


S?  LET  vr*=FACT  *SUHV<2>/CCF/LAP6F'A«*?-rTX*»2 

5 7  IF  VT*  If  \J 

K«  PRIM  1  LIKE  .ITH  VTX 

CS  THUS 

Tnf  UnLt  !'  .FI.TRl’TC.  VTX  =  . 


-YSTT^b,  °rLE«SE  ? 

DEC  *  9f 


O 

•X 


r 


a  r-4  K>  (\,C.(\j 

c-iii  <r  in  m  in 


cm  a:  r-  r-*o 

c  <v  m  «■  fo 


OK.CC^ 

4  «-K>in 


*-h0'  c.^n *«m «-><  r-cr .-hit  x  ccc  cr  o  oin*o  <\jk> 

cjcg  «r  1/  »om  ^ «■  cvforocvjKiu  *•■<«-  x*  in  c\.cm 


-*x*-»(T‘4n»-'Co<^  h*r*-KU)*~»cs<h*r^r»ricr'a(GO^'«vj*-ir-*r 

'"'JCNi  lOK,  f\  IT  «-i  u  riK.FOlfl  <\J 


cxarir  a:  a:  cn 

LuLiU'UIL-U.UJUJLJU-iUJ  W  UJ  LJ  L*J  UJ  UJ  UJ  Ui  l±J  UJ  LJ  LJ  Ui  O 
JJJJJJJJOOC  JJOCJJJJJJJOJJ 
ffiXJU.aiQLCCa^C^UjL^Uj  ££  £DLJLJCC£Ij£LCDODOQCOLja3CO 
D3HK3DDDDDDHDD 
OOCC  CC'C0^22  C0220000000I-00 

Q'-JC  ODC-OUmhh  OOHHOOODCOOHDC: 


-IV  CvU-tf  ^  •— »C  'K.  <Vo<*lT^U 

Cj  KUV  K'  C* 


i  *  <*  x  r-  * . 

p*,C>JK> 


.  J  Uu  COQC  QGul-i  O  Cj'-jC  O 

•  a  uq;  •  bora'a  «  c*  a.  a  a  •  oactoc  a.® 

r>u  foe  c  err  o  cocceo  o ooc  oe 

*;  j  SJ2  jjja  233332  ajja 


U  U>WLl 
-JO_J  -J 

a  a.a  aj 


a.  a  u  Li  u  o  a 

■-4  4  4  4  4  <J  4 

>  >>  >>>> 


U>U<U  U»  U.IL.UL1  Ui 

.JO- J-i  J  J  J-J  -j 

a  c  a.  a  ce  a  a  .’  a  a_ 

4  4  4  4  4  4  <  *J  4 

a  a  a  a  a  a  a  a  ct 

«<l<4  4444  4 

>>>>  >>>>  > 


i.  >  ..  U  U.U.UJO  k'luU.U  UJU.U'L.1  14 

►“  >  >  >  K-  >>>>  »->>>>*-  >  >  -»  >  >  H 

*r  — *u  ^ 2 uj»-* *“•*-»*-, Uj»-* 2 

u  /I*  oiL/lUJ*:  W0iy0i/)UJ  14  '/)  </>  */l  GO  u .  •  i/>  i/)  1/7  </)  £  l/>  Oi 

x  ur*"u.Li:  *-o.aaa  Jl  ixxaa  X  *~«a.  u.  x.  a  *-»a:  2 
->  O  ►  33J>-33D3UJ033PD>-D333KD3 
C  »_  e  JttLUlCOUOUUI.  ibOUuDLO 

l  uj'  4,1.1  j  C  uiUUuJ<^UJUl>JUiA^UJuJUU(.  Ub. 

-  j  .*  u  4  «  uicllOi'i—  *iaiiait«iau  aira.aa*' 


.*  U;  JU</)*  • 
■  V  0  V*  -la:  ►  • 


Jl  U. 

<x  aoi  •  >- 

►-  u.  »  ;  ►  xc.  u 

L4.iQn.4j_j*-  x  -x 

u  OUl  u.UILiU.*' 


T,  *►-:>*  »OX 

&  J.  *-  lI  >  a  x  4 

uodu  Huimi 

2*  a  (/,^>-r>jxx 


B-17 


TNT  7  L  INr  S  WITH  K  ,  f;  ,  ANS  WE  R  ,0  I S  T  .T  vpe  <  PI  ,  P2  ,  am  ,  $  TD  V 


DI^TPT8UTT  C  ►'  FUNCTION  rCR  THU  *•  TU  LARGEST  rf  ..  .........  .........  p  v«s 


CNi  OJ 

♦  *  X 

*  ♦  o 

X  l/) 

U  Ui  » 

I  I i 

c 

•  «  ■ 

K)  K)  »  *  • 

V.  VX  «  < 

M  PNJUJ  *  ' 


O  O  *  «  ■ 


f  fi  r\j  -*■  qhuzm 

|  u.  ~  II  II  0*1  •o»-»  O  *CJ_J 

|  II  KjCVJU.  a.  *JZ«Q  «  «<M 

|  c  «  UJ  »U  Ui  ?  Q.  t3  •aCMSI 

If  C  “5  OUJ^C  Cj  D£  II  II  2  *  **  •*  Oh 

—  -)  QwqujO  U.O<-*hU.O<hi.— 
j  •  •  •  *4>  1  *  L_  V)  N.y^WO  *  P  *  U_ «;  </) 

|»rtr  HO  XOC  —  ||  IVNKHWWH  3 

*  ii  n  ii  iiojujO  houjoiujoo^s  >>u  n  >>  n  uai 
III  ^#.«^|-U)>KNCj«J>*COU.U.aOli.U.*JOQ.h 
N*-*cmk>^  o  i  u  <  a. 

tjww<»v<Hh«  MH  X  JKHhhKHHHa 

J>>>»  II  Ui*  II  UiU  »-  .JUiUlUJUJLJUiUUdOU. 
fciU.U_U.U_  *D  -I  *-Jm  o  <JJJJJJJJU-H 

I  a. 

l-l-KI-l-OC  OS 
LuU^UiUUiO  OO 
UJJJJ^  -tli 


OUiUi  UJUJUJ  UJ  UiUJ  o  u. 


src*  f/y 

0*H»  3 

O'  -•  X 

h- 

«  *11 

*  c.’iouMi/i 

♦  u. 

*  Aa  Maz 

•  O  003UJH 

«>Zu.co:>J 

«<£  Ui  o 

a  X*  <Ni 

JU.  fr- » 

<*M  O  H 

*  cl  r 

*  Oh 

*  ooc 

*  JCL 


HH  Nj  ♦  fsi «  XX 

S  U_  JOJoI  I- fr¬ 
ee)  Ui  *Ui  «fr- 
Ui  o  OCMOIVM  zz 
m2  •  n*  *  w —i*-* 

IlH  U_  ‘H'U.H.U 

CJ  -J  O  *h  *»r  #v,  > 

<  ^fr-^fr-ui  ore 

-J—  i-  >ac>o:^  :d 
U_  Z  XOXC5*-i  XX 
•fr-  M  3V)3MJ  DM 
UZlAC  O  II  O  II  XX 
ZmDO.  IIK|Hx«- 
3aeX  vjnqxq  xx 

o;a.fr“Z  >~uic>uiv>»-  «* 

*-  0«*  ZlAtU. 

u  fr-ouJUiuiwSxu-0 


I 

ir  .io^rvnvDr^  orrr  <  «-«<v»0  4nr.x>r-  ecu'  c'»r*<<\,K'  ♦  it  sTf^ao  ♦  on  o *r  ♦  ip 
iPiniTiinipiru  iT)uOvovO'OU)Nij^5sO'X/«£fr^^f,-r^H>^-h*rs‘^ *  coctcoccco® 

I  • 


vf  f*-COOCO'<^’**CVK>«-lf  Q 
00  00  00  ui  00  O'  O'  ON  ON  ON  ON  z  2 
0.  *■  < 
G.  Uifr— 

3  XV) 


M 


“a'n  poi'T’nt  CficTrTwcrPTpTri#*.  roRppi**^svsTrMs»  p^iEtsr  2.1  *aGf 

OPTIONS  =  SEQUENCE • IO*r UL^wk #vRfF  ,N0£  vPLT^T  «T?ti*r?  r*  DEC  19*4  "c:k0:?3 


♦  f-  r~ 
ff4  *fl 


DC-  (1  H 

2*  *-  cc  ^ 


iT^ct  \flh-r- 


O  vfl 

ir>  u>  «  ♦ 


in  in  «  ••  < «  K>K>h»o  HOiO  or-  h-r*-K>  m  ♦  mo  <c  c  cr  r-o 

rH  r-  cmdct  cr  mr- f- r-  -ujinin  sc<-  cgr-  ♦  «  r-  rm  in  m#-  ®f- 


♦  ♦  ♦  *  * 

+0'tf  r^ir  cr c. coir k>%dsc^cm«-»  «h  »o*- <j *Hr->  (r»o«-»*-c\j  mr- 

^^r-ia^cr  r-vflvo^-annm  era:  oj<rr-c\jK«*  ♦  r-crcr O' ♦  irc\jfv <r  ♦  s £ 


F'*(TK)s0C\lO(\jC'<Sh>(00'h-<0(h^lfiM/^^M0'^K)(Mti)HK)HinHMl0^l0inC0MU)Hf>4t0MH 
rH  *fi<r  ifirJU'O1*  vD  vUsL  VL  -h  KCN.UC*-*  ♦  vD  MOMONMTO'^^rt  O'  •-*  r-Uf  <T  <•  \i 


UiUUU  UIUJ  UJUJUJ 

»-  ff>  OD  CL  CL  K  CL  CD  03  02  CD 

yOOOOUOOOOO 

hCLQCKCOOOD 


actccc 

It  I.  IL.II .  II.  It. 1 1.  I 

JJJOUOJ 

(ECDCDUJUiUi<X3 

333Khl-D 

0002220 

QQCjmmhQ 


ac  cox. 

UJ  UJUJUJ  UJUJUJ 
JOJOOJJ 

cd  uj  oq  uj  ui  oo  ao 
Ol-OI-H-OO 
0202200 
o  mOhhOO 


OCX  OL 

Uj  UliiUJUU  UJ  UJ  UJ  UJ  UJ  Ui  Lu 
«J  JJJJJ  -J  OO  JO  J  J 
CD  0DflDfO(DQ3h-IO  LUJ^UJlSCD 
3  33333X3  »-»-3»-33 

O  00000^0  220200 

O  OOQOChO  hhOnOo 


<£!(  C. 

vcc^«t  a) 

r- 

CV  ^  Cl  rH CVK 

♦  r* 

40  uu 

r 

a  sccv* 

m 

o- 

h  ♦K. 

0  CC 

< 

m  in  in 

<r  -rot 

IO 

♦  OJ 

0 

■a  IT*  in 

m 

o*m 

0  0 

CJ  000 

0 

000000 

OO 

aoo 

c 

OOO 

0 

0 

00 

a:  ua: 

aca.Qc.ocac 

DC 

CL  CL  CL  CL  CL  a 

a:ac 

oca:  cl 

oc 

acacoc 

a: 

oc 

aco)oc 

oca: 

OftC 

c  aooo 

0 

000000 

oc 

000 

c. 

COO 

0 

c. 

a:>-a: 

oc 

i«o 

303J 

3 

33313133 

33 

333 

3 

333 

3 

3 

«<(/)« 

33 

UJ  Uj 

U  UJL./UJ 

UJ 

U  UJWUJUJU 

Uj  uj 

UjUjL 

Uj 

UjUJUj 

UJ 

UJ 

UJ  Uj 

J  J 

J  JJJ 

J 

1  1  »- i_J  < 

-J-i 

JJJ 

-J 

— i— 1 

-J 

-1 

-J.J 

a  a 

an  ma  oj 

CL 

(DOOftfiLa  CL 

rap 

a.  a  a. 

a 

a-axi 

a 

ff 

a 

or 

«a 

♦ 

♦ 

u  *— u. 

aJ  K' 

►  •  JHHH 

•M 

►•HMHHH 

►* 

H 

*-H 

-JKJ 

hh 

uma 

a:  cd  a:  acoc 

a. 

a:  a:  a:  a.  a:  a: 

oca: 

a:  oca: 

CL 

ockol 

CL 

CL 

aioaj 

cl  a. 

<«<< 

« 

<</>♦ 

>*-*> 

>•-«>>> 

> 

>>>>>> 

>> 

>>> 

> 

>>> 

> 

> 

*-*CL*-« 

>» 

a 

a. 

QC301 

UJ  <  Uj 

uj  «l  L  UJ  UJ 

UJ 

UJUJUJUiUJUJ 

UiUJ 

UJUJUJ 

UJ 

UjUJUJ 

UJ 

u 

<co< 

UiUJ 

>» 

>»» 

> 

>»>» 

» 

>>> 

> 

»» 

> 

> 

>  > 

>> 

h  »^uj*«  hhhLjh  Ui  ►*»-*  UJ  *■*►■*  *-• 

«/>  ~j  in  z  co^mco  </)£</)  (/x/x/nn  </></)«:  2  wvizcowto 

a«ut<H(L<iaaQ:^«  KttKKKKH  jjhmcaHaaa: 
3a-D»“3cij333h3  333333HLiUl-33h333 

UOUOUOLiUOOU  UUUt)UU33}lS3U03UUU 
UJ  JUJC  UJ  JUJUIUJOUJ  UJUJUJUIUJUJO****  O UJUi OUJUJUJ 

Q- oa.ua  oko-kko:  oca: a:a: a: ck: a;  jj a: oca: a. a: oca: 


UJ  UJ  ►-*  O  U*"H 
</></></>  Z</>Z</>  JUJ2U)W 
(Lococi-iflc^a.  <M«Ma:a: 
D33KDH3  DJS.H33 
U003030  oaooou 

UJUJUJOUJOUJ  JljOUU 

GcocccxococQc  OHoaa:  a. 


U  U  CJ  </> 
l  t  H  u, 


♦  ♦  k 

• 

UJ 

>- 

« 

0 

UJ 

UI 

a. 

J 

IT 

u 

•J 

z 

u  > 

►* 

>- 

WU- 

a. 

u 

• 

U. 

3 

Oui 

u 

CL  IaJ 

►- 

>0 

a 

O 

• 

U  CL 

• 

OC  u 

n  U3f»- 

l 

ua  -J 

• 

00 

0 

■ 

UJ 

u 

•  H- 

03 

k2 

2  -ro 

< 

3 

— i  • 

00 

Z 

• 

2  • 

»-></)> 

2  • 

•3 

<  c 

2' 

ouu-i  -J«/i 

cr 

•  j 

u 

0 

♦  c. 

xNarOffii 

3a  o.  •— a 

Uj  ►* 

2QOOUI 

H 

xrsjl3> 

220 

0‘ 

■— 0 

UIO 

rH  CW 

0001-33 

at*-  ►* 

u.u» 

X.  0 

iuuuaoyyu.L  u 

coiH3^kJJJiir2aa  a 

<3  00</)C/H^</)10 

►-003X1 

B-21 


V. 

u 

a 

o 

CM  ♦ 

2 

a  a 

a 

oc 

a 

*) 

a 

a 

• 

•  , 

a 

a  nocj  O' 

c 

a  a  CM  CM  a 

CM 

a 

a 

o 

i 

V) 

a 

a 

O'  O  C3  OaCMLTIP  VT  Okf) 

m 

O' 

X 

=> 

< / 

2 

«j  o 

DJUJ 

a 

-JO 

2 

K)coifiaaHOOHHinCa 

aaCMa  aa  aa 

a  u 

a 

t/> 

• 

3 

K  fV 

oo  & 

a  act  a 

Vi 

Q  r> 

L.  Jl.  il.  Jl.  JI.H.  11.  1  I.  U.IIaI 

K- 

a2 

ao 

a 

JJJJOJCOOJJJJJ 

o< 

h- 

o 

fli  ao  to  cd a  cd  a  aa  oo  <s  oo  ®  cd 

>- 

or.  X 

aa 

o 

DDDDH3HHH3DDDD 

<1.  K 

&  X 

* 

ooccroz 2200000 

-J 

«a  b- 

QOODhOhhmooudO 

Ul 

ir. 

3 

a  </j 

«t  ct 

►-a 

C  i 

or 

a 

u 

a 

* 

X 

v> 

H-  V. 

a 

L-*“ 

*-*  — i 

c 

X 

-*  a  • 

H-  » 

<j  *•* 

*  U  —  CM  aC  H  ^  h  U  U 

a  <_>c« 

u> 

H  N 

a  => 

• 

4l  w- 

uoz 

t-. 

a  u 

•  (/»  » 

QO  OO  CO 

i~ 

aa  •  *aa  •  *aa  • 

a 

Zj*.^ 

oo  c  oooccoco 

K 

►  O 

a  ? 

33  2  23322332 

»-o 

-  * 

a 

Sc  O  u. 

a 

-IX 

(J*-  X 

o 

L- 

»-  <i— 

l_l>-  w 

t/jr*. 

a*-< 

hL’h 

S' 

l  lj  a  a  aa 

3  MJ 

u7ir 

aa  aa  a^ 

at 

GD 

a 

c-cc  a.  CD  COD 

LI 

I-a  a 

««  «*.  X  «  X. 

a 

a  *-* 

a  • 

OJ 

a 

aa  aa  aa 

►  *- 

a«< 

0,44. 

a 

a  a.  aa  aa 

o  • 

a 

a*-*  l 

a 

CM 

a 

*«<•  -*.< 

«  K- 

a 

OX  1 

03 

X) 

>>  »  >> 

o</» 

>  • 

o#-  ~ 

«* 

2 

a 

a 

a  o'; 

o 

a 

» 

aa  aa  aa 

-J 

• 

►-a 

a  ►« 

a: 

a 

a 

a 

>>  a  H»HH»H 

a 

OJ 

<a 

■a  »a 

< 

3 

a 

wauiU2UZNN22NH2 

x 

a 

-JOJ 

it  nD 

> 

2 

a 

► 

oooo  Z2uzuv)v)ua^v>u 

D  -C 

U>  LLi 

• 

a 

qckhhj  axaaxxaai 

* 

*  *-* 

a 

a 

» 

DDHHDH33DSDD33 

2 

<r~* 

oa 

UZ(£ 

a 

X 

U033030UUUCUUO 

• 

a 

o< 

iuih 

o 

00 

aaooa oaaaaaaaa 

> 

%  iO 

a  ~ 

a 

aaaa<ta<aa«4»aa*i 

u 

• 

a 

L 

a 

V) 

V) 

i* 

li 

x  oo 

2  » 

X  Hi 

X 

*-<_» 

a  or  1 

M  CM 

3  >0 

C1 

V 

• 

c- 

«ou. 

a 

•  a  • 

* 

00  A* 

oo 

oo  *■*  * 

Cj  sl. 

(X 

X 

a* 

JIW 

« a  m 

a  I 

o 

e-.t 

aa  » 

♦  a 

X  xo 

o 

*4 

<M*  «- 

a  i  « 

u 

w. 

-Jt> 

UliJ 

13  *  u_ 

»-■  ON 

„ • 

*— 4 

^  o 

r 

2  2  O 

.B  • 

• 

O 

oo 

7f 

•  V1C  mUCJ 

axa 

J 

• 

ao 

C  2m 

a  «  *  ♦  •  «z  Vh  00 

H 

2 

<  4f 

a 

OCJCM  kZO.KL.IHH 

• 

u. 

o*-> 

a 

2  *1  •'  ID 

UJ 

c 

a 

UM 

•  Ht.K2L.UHNH  • 

CJ 

* 

22 

2  ao 

2  33  11  II  II  JU3L.Z2 

2 

Oo 

■D  C t 

•Z741K 

aaxoca 

a 

2* 

a  c  u 

3C  X 

a  oo 

a 

C' 

N- 

aaaaa 

IHHK 

L5 

► 

a*- 

Vi  •*  u 

aaaaaaa 

aaa* 

L. 

¥- 

2 

►-•  o 

2aaaa_i«-» 

Oja* 

VJ 

l. 

JU 

+  r  iX 

M 

•■ 

a  a 

*-  *•  c. 

u 

c. 

II 

J* 

1  •  •  • 

•  •  •  • 

a 

2 

a 

U-» 

•  •  •  • 

•  ♦  •  • 

c. 

a 

00 

00 

a  *-* 

7- 

a 

►-o  o 

<  5 

aCMK;  *  U 

v»  a  a  CT 

«  •  ^4 CM HMTtf  \CHCCffOH 

2. 

V>«J  • 

nnhNM 

« 

a  •  2 

2 

02  a  -*<v 

•  3  o  3Dn{\j  a 

’ 

«Sau.kKZZZQ.(LKKN 

I-Xj 

O' 


(U 

JUJ 

Jl~ 

:u> 

» 

Cd 

Sr 

m 

« 

♦ 

<M 

H 

i 

5* 

CM 

*Si 

* 

« 

%» 

«■» 

k 

H 

• 

♦ 

— i 

Jr 

«i 

• 

Uj 

cm 

0 C 

«* 

V 

u 

V* 

• 

Jl 

* 

_J 

1 

U 

«a 

^■4 

U 

U 

♦ 

u 

* 

u 

v> 

CM 

Zj 

H 

• 

1/t 

• 

*■* 

CM 

* 

I 

3 

CM 

t 

Jr 

* r 

3 

>-< 

• 

ar 

_J 

CM 

u 

« 

a 

♦ 

• 

* 

» 

cm 

> 

H 

u . 

ZJ 

►H 

WO 

C- 

s 

a 

o 

* 

+* 

*» 

o 

» 

u 

u. 

♦ 

U. 

-X  • 

H  • 

H 

U  1 

•  -J 

«  -J 

* 

fr 

«  « 

H 

** 

cm  us; 

CM  UJ 

</)  V> 

• 

•  *3 

•  a: 

*  N-» 

* 

*-*  O  «  CO 

CM  O  « 

a  o 

X 

3  OS 

3  OQ 

•  ©  • 

t 

•  a.  c< 

X  • 

U  1  u 

u 

<i  —  0  _jhh  *  •**  *■*  ♦  -JCM* 

•  4 

; 

U  x  M 

«m  u  c».  sr  • 

i-  ar 

c 

_j  •  _j  as 

u  cm  j  •  _j  uaaic 

*'  •  •  >  «  O  OL  ft  Uj  >  '  C  tlUiV.H»~* 

vJ 

OwOir>^-«,-»  OUJHHH>UlOUJO.#i  r>»-*  OUIHU>C*MH 

ahw  •r»ii  iiciu»-m 

*M  *OkjH(0  O 

• 

Mt.JI  *  II  «J  •— ~>“Ci 

•  «r*jCMII  J  ^t-il~»  CM  Z 

Ul 

»-«<  jn  z>  ft  j  caiiuo*  M(  an  ii  -j  c  i  n  uj* 

iizo 

O 

00*(BO^I-X  JO*Xu« 

ooatuoux  jooc sc  «j 

fitt  u 

XT 

*  u<  h  sc  u 

a.*  Uj  h  Jr  UJ  Q. 

3 

'*> 

x  **  h  hh  x a 

O  XH  HU*  ZK  OHHU 

J) 

U  *»  Ul  C  Uj  W  UJ  U  HO 

OU  HUJOUJUU  HO  OkJUCi 

J 

— -  u  _J  ^  _l  -J-J^  Ok. 

_*•-  OJCJJI-  OU 

Uj 

m  m 

•  • 

•  »  • 

Vi 

*  ^ 

0  f  • 

♦  »« 

#-«  »-• 

•  ( 

c  *•  n 

II 

-J  -i 

-I  -J 

-j-iar 

•  • 

•  • 

•  •  UJ 

IV> 

»o  «  <r  -  •  f  *  *  w  -cm*’*#  i/WHcpo*c-*«-<cs.»^ft  tri'XJH® 

w  1/  ui  x  ^  M  u)',“Hh^i“-hh^  hh  cr  «  cseDffgoss: 


!0«JS  -  SEQUENCE. ID.SUECHK. XREF.NOEXPLI ST.TFtCE? 


O  CL  ¥  i 

♦O  OJtfl 


(T  mO'-CT'  0  O' 

— »IPCM^  in  ♦ 


\D* OCMCL  vQ  00 

♦  ♦ 


•  * 

mo~  vxjcj  o  •-«  «<ra'  cm  oj 

♦  ♦  <r  u  ip  nl 


CM  CMCC  iff-  «- 

UP  •-*♦*-  K)  CM  ♦ 


«  ♦  p-  a'aa' 

♦mm  ♦  mm 


lj  vO  m  »-*0'oc<rosjom»ocr*-‘a'#' 

cu  if-  m  K> ♦cyfO<r  ♦  ♦*£  ♦  k) 


C0£T>  nO'K)lf)H  4)  a) 

c\  cm  ifirju-ir^  ♦mm 


z  m«~4«^K)(-4  0\aoioc.  voa3fom«Ho  •-•(r'inv-tfooaxPloo  aocM^O'a>%ciom 

tf  <**|fl|n  <V  ♦  «-«<MP-  fOlOlP  ♦  CM  If-  ♦  CM  <r  i  v  mu  lu  ♦♦mu 


LJ  UJ  UJLJLJ  Li  LJ 
Lw  CJJJJJJ 

a  uiajaimaLmcD 

O  K3333D3 
»  2000000 
*-000000 


oc  a;  ac 

U.  UJ  UJ  L- DLi  DbJL.  ijGUjLJLj  L.LJ  LJ 

O  O 

UJ  UJ  QOCLflQUJCDCLCOCQGOCDSDCOIACOQLi 

►“  33DKDDDDD33333D 

2  2  COOZO OOOOOOOOO G 

•-«  QOOHQOQLOOCOOOO 


Ki 

a  cm 

r 

K 

HIOh 

♦  UJsor'^ 

c  Q  cm 

~<<M 

rO 

CM—CMinrO 

CMCJK) 

Lj 

CO 

O 

00000 

000 

QC 

oca 

• 

• 

•  mac 

a  a cocK.ec 

ficacac 

G 

oc 

0 

O 

o>-o 

cooco 

OOG 

31 

39 

z 

2 

zwa 

73333 

333 

Ui 

L> 

L_  U' 

H-L) 

U'UJ  UjLjLJ 

U  ULi 

-J 

-l-J 

3J 

rJjJJJ 

OOO 

U 

cr  qq 

OLCL 

0  CLCLfflQD 

O  CLOU 

♦  <♦ 

•** 

►HM* 

QC*“< 

HmH 

cr 

aa. 

►-QC 

QC  QC  Q.  LC  OC 

a  QC  CL 

< 

♦  ♦ 

*-♦ 

> 

>> 

<> 

>>»> 

>>> 

L 

UJU 

KUJ 

U-LJUJU.UJ 

UJLJLJ 

> 

>> 

I- 

►~z> 

»>» 

»» 

Q.  UJ*-UjUJUJ>-»-  Z  Z 

>-  i  M2  22MV)  UJ  W 

*«qc  oc-j«j_u«i-iaL  2 

>-3»-»-V-33LJLJLJLJLJLJ3  3 
^>U  333OU/CD03C0  CDCDiOO  O 
C  ijOOQUJUJ««  ft. 

a  acLoctao:  jjjjjJc  <i 


Z UJ^UJUiLJ »~*i  ihiHMU»  M  <m 

ui  z  co  z  z  z  uo  </>  </>  to  co  z  cn  </>  </> 
*.  «  a  •-«  m«  im  ac  cl  ac  a:  a  ►-  ac  a:  a 


OKU33DULUUU3UUU 
CCUIUJO  OOUJujUjUJUJOuJUjUJ 
*  ft-  ac  ac  oc  a;  a.  a:  qc  a  a:  ft  <r  a:  ac 


2  CO 

U  L  .CMCLfOW  K- 
3^KXX>XQ 
l|  II  II  II  .J  • 
LJ  3  K)*-J<Vr-  U. 

a  cr  x  x  *t  >-  c 

U  K  u 

O 

►  LiLJLJLJC  • 

G  jJJJC* 


r-  o,(T  r 

if- mu  %fi  \f  *0  «ii  .fl 


c/.  a:  a  c  c.  .-1  cm 

CD  or  ♦  iPt/iN  aj  3  3 

♦  DUIU.U.U.  I  JJJJJJZ  Z 


•  <C2  O 

mi  LJO«  HCj(OZHHrg 

O.Q.craMhDixxxxxx^v 


B-27 


-  C  U  T I K  F  FOR  T^E  STANDARD  NORMAL  PROpfFTLlTv  INVERSE  cU».CTICw.  POUTIf,E 

TURNS  ThE  CUA\TlLr  rF  ThE  NORMAL  C.C.F.  GIVEN  THE  CUMULATIVE  PROBABILITY 
'^“.PROP*.  r’rtrERE\CE:  A**S  ♦ H £  I rc ^ 0 F  CF  MATHEMAT  ICiL  F  l»CT  I  OHc*  . 

'  :T.  euREAU  OF  STANDARDS*  (F.  F33>. 


•  V* 

C\  J. 

O' 


_J  Cl 
UJ 

a  u 


co 

c. 


2  • 

2 


C  Ci 
C 

L 


•—  a, 
_  c 

►-  a 


c  <1 

U  -L 
u.  'A 


_>»- 
1  ^ 

Do 


*«  & 
—  O 
a:  co 
O 

*-  o 


Y~  «*  X 
2  » 

20  -  ■  *- 


»-  O 

• 

X 

h-  li 


-J  Wj 


a: 

CJ 


ft:  r-  Ch-ir.  a»cr»fl 

u. 

o 

co 

qc  *  ♦  «  * 

UJ  stJ  O'-  tf,  <t  >i,NN 

CO  ^  r-« 


c\4 *-ir- 


a  a: 

Cl  Cl  C  C  Ci  Cl  Ci  Cl  UJ 
UJ  JJJJJL30-IJ 
O  COCCffJCCCCCuJCDCX 
O  D3Z>3  3l>-l-=>3 

s.  ooccozzoo 

QOOOOmhUq 


u 

Cl 


Oj  <r  —  cv  cr 


a  oao  u 

ac  ft  ex  q:  «a:  • 

c  ooocoo 

3  33132312 


Cl  O  OftJ  c 

-I  -J  o  _J  _J 

an  er  a  od  ex 

«*  «  <a  <  < 

u.  ft  a.  a.  a 

«  <<<  «a 

>  >>>  > 

W  cj  c  cj  Ci 

UJ  >  >»>►-»»- 

a  huUj^hmZhz 

>-  C0  2i!CO<OcOClCOftj 
►-  oci-H*«ftcc  iriai 

UD3UUUOOC 
uocuuuKuja: 
oca  a  ftft 


ii 


' J  •'  I  •  *-•  w  OjCo  HK 

—  *—  *1  3  <  i.  <ftC  CO 

i:.  >*  hqq  m<xij  X*-h 

I  tC  2c  2  ft  •  *  —  ft  ♦  l-c 

O  <  •  »  dUC  H^  (  QKitHH  « 

00  L..X  .3  CO  M-  |i  II  z>-  II  N*  30 

2  X  HMV.I01  OOOft  2 

2<a  -x  x*-c  >«ftOo-rQCftc*2< 

O  ►-  t.J  >l  ft  ftU.  <COK-lu  <  OKJ 

—  <r  «  L.  OSK 

—  *—  i  l-UJ 

C  T.  JU  u  K-UJUJCI  n  u/co*  Uj* 

*  ~  2^-  OOOo*-*_l-J«  ft • 

:  >  -  •  u 

*  *  c  C3  c. 


*r  r  -i'.t*.  a  c-r  ^<vio  «  u  i>f  f-ap'  o 

— <  • - ^  '  Cy 


O  ft 
2c  C 

<  IH*- 

o-jo  ao  x 

ftft  •*  X  ««  3 

«1  Cft  CJC  ^2  to  X 


B-29 


i  3a  JS*+ 


mm* 


-.19 


P  F  L  C  «  S  C  2  . 

c  ore  i9c« 


*iCM^4T  f*  a?CTC  *-*C\|K)  ♦  tfi  \£  ct  Q-  .  ■>»-  -•  Cj  rO  ♦  ID  s£  r-  OO  CT  C  »Hf\lK)CinuJh  |3*00HQ: 

<NiC\.<\CVCy*OK.»OfOK»K>K>»OK)fOtf  «-  <9  <t  ^  *3  <t  >t  <r  Kf)ir)C 


^O  T'f.r  "  “  F  L  C  T  £  .C-  A  CAC’  f!,'?''RIpT  M.K  o  p  i  u  r  <;YSTE£T,  BfLEASE  2 P£Gr 

CATIONS  -  r  E  t.LF  NCE  .  TO  Chi  .>  REE  .\G£  xPL  lrT  .  Tf  £  Cr  ’  ‘  c  DEC  1  96  *  '  °  i :  2  ,  :  ’  3 


END 

FILMED 

3-85 


DTIC 


